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GEOLOGICAL-SEISMOLOGICAL  EVALUATION  OF  EARTHQUAKE  HAZARDS 
FOR  APPURTENANT  STRUCTURES  AT  GATHRIGHT  DAM,  VIRGINIA 

PART  I :  INTRODUCTION 

Purpose  and  Scope 

1.  The  purpose  of  this  investigation  is  to  define  the  maximum  potential 
for  earthquakes  at  Gathright  Dam,  Virginia,  and  to  provide  time  histories  for 
earthquake  motions  that  represent  the  cyclic  shaking  that  would  be  felt  in  the 
free  field  on  bedrock  at  the  damsite.  Ground  motions  defined  by  this  study 
are  for  use  in  the  engineering- seismic  evaluation  of  appurtenant  structures  at 
Gathright  Dam. 

2.  The  investigation  includes  both  a  geological  and  seismological 
analysis  and  consists  of  the  following  parts;  (a)  an  examination  of  the  local 
and  regional  geology  with  an  evaluation  of  faulting,  (b)  a  review  of  the 
historical  seismicity  for  the  area  under  study,  and  (c)  the  determination  of 
the  maximum  earthquake (s)  that  will  affect  Gathright  Dam  together  with  the 
attenuated  peak  ground  motions  at  the  damsite.  The  maximum  earthquake  ground 
motions  specified  are  in  accordance  with  the  requirements  mandated  by  ER  1110- 
2-1806  of  16  May  1983. 


Study  Area 

3.  The  area  covered  by  this  study  includes  that  portion  of  the 
southeastern  United  States  in  which  earthquake  activity  has  the  potential  to 
affect  Gathright  Dam.  Gathright  Dam  is  located  in  Virginia,  near  the  border 
between  Virginia  and  West  Virginia  (see  Figure  1).  The  study  area  is  in 
general  limited  to  the  region  contained  within  a  circle,  with  the  reservoir 
formed  by  Gathright  Dam  at  its  center,  and  having  a  radius  of  approximately 
100  miles  (160  km).  The  study  area  includes  portions  of  Virginia  and  West 
Virginia  and  incorporates  Giles  County,  Virginia.  Giles  County  is  the  site  of 
the  second  largest  historic  earthquake  in  the  southeastern  United  States.  The 
Giles  County  earthquake  occurred  in  1897  and  was  felt  over  much  of  the 
southeastern  and  eastern  parts  of  the  United  States. 
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Figure  1.  Map  showing  the  location  of  Gathright  Dam  and  physiographic 
subdivisions  in  Virginia  and  West  Virginia 
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4.  Gathright  Dam  is  located  on  the  Jackson  River,  approximately  10 
miles  (16  km)  upstream  from  Covington,  Virginia.  The  dam  is  located  in 
northern  Alleghany  County  and  the  majority  of  the  reservoir.  Lake  Moomaw,  is 
located  in  southern  Bath  County.  Lake  Moomaw  is  approximately  12  miles  (19 
km)  long  and  ranges  from  less  than  1/4  to  1-1/2  miles  (1/2  to  2-1/2  km)  in 
width.  Gathright  Dam  is  a  rolled  rookfill  embankment  with  an  impervious 
compacted  earth  core  and  a  concrete  cut-off  wall  in  the  left  abutment  for 
seepage  control  (U.S.  Army  Corps  of  Engineers,  1983a).  Construction  of 
Gathright  Dam  was  begun  in  1967  and  was  completed  in  1980.  Gathright  dam  is  a 
multipurpose  dam,  providing  flood  control  and  recreation.  The  dam  and 
reservoir  are  operated  by  the  U.S.  Army  Corps  of  Engineers,  Norfolk  District. 
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PART  II:  GEOLOGY 


Geologic  Setting 

5.  The  study  area  is  situated  in  the  Central  Appalachian  Mountains, 
approximately  at  the  junction  between  the  southern  and  central  segments  of  the 
Appalachian  chain.  The  Appalachian  Mountains  are  a  major  physiographic 
feature  that  extend  through  the  eastern  United  States.  They  begin  in  Western 
Alabama  and  continue  into  the  eastern  provinces  of  Canada.  The  Appalachian 
Mountains  are  dominated  by  intense  folding,  with  numerous  faults,  and  are 
composed  of  a  vast  variety  of  sedimentary,  metamorphic,  and  igneous  rocks. 

6.  The  Appalachian  Mountains  are  subdivided  into  geologic  provinces 
based  on  similar  rock  types  and  stratigraphy,  structural  features,  geologic 
history,  and  similar  topography.  The  major  subdivisions  in  the  Central 
Appalachians  are  the  Appalachian  Plateau,  the  Valley  and  Ridge,  the  Blue 
Ridge,  and  the  Piedmont  Provinces  (see  Figure  1).  Gathright  Dam  is  located  in 
the  Valley  and  Ridge  Province,  near  the  Appalachian  Plateau  boundary. 

7.  The  Valley  and  Ridge  Province  is  composed  of  intensely  folded  and 
faulted,  mostly  unmetamorphosed,  Paleozoic  (600  to  250  million  years  (m.y.)) 
age  sedimentary  rocks.  In  contrast,  the  Appalachian  Plateau  is  composed  of 
relatively  undeformed  Paleozoic  sedimentary  rocks.  Crystalline  rocks  occur  at 
the  surface  to  the  east,  in  the  Blue  Ridge  and  Piedmont  provinces.  The  Blue 
Ridge  Province  is  composed  of  thrust-faulted  Precambrian  (before  600  m.y.) 
basement  (igneous  and  metamorphic)  rocks.  The  Piedmont  Province  is  composed 
of  highly  metamorphosed  and  folded  Paleozoic  and  possibly  Precambrian 
sedimentary  and  igneous  rocks . 


Tectonic  History 

8.  A  brief  summary  of  the  tectonic  history  of  the  southeastern  United 
States  is  presented  below  as  it  aids  in  understanding  the  present  geology  and 
seismicity  of  the  study  area.  The  geology  and  structure  of  the  different 
geologic  provinces  identifies  a  complex  tectonic  history  that  involve  multiple 
periods  of  deformation  during  the  past  600  million  years.  The  geologic 
history  includes  wide-spread  volcanism,  metamorphism ,  and  several  collisions 
of  the  eastern  North  American  continent  with  other  crustal  fragments  (island 
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arcs)  and  the  African  continent  (Cook  and  others,  1979,  1981,  and  1982, 
Hatcher,  19/1  and  1978;  Rankin,  1975;  Van  Der  Voo,  1979;  and  Williams  and 
Hatcher  j.982).  The  major  tectonic  features  identified  in  Figure  2  were 
produced  during  the  Paleozoic  as  a  result  of  a  series  of  continental 
collisions.  Associated  faults  shown  in  Figure  2  are  from  Price  (1986),  Calver 
(1963),  Butts  (1933),  Lowery  and  others  (1971),  and  Johnson  (1977). 

9.  Williams  and  Hatcher  (1982)  proposed  a  model  for  the  Appalachians  as 
a  mosaic  of  tectonic  terrains  or  suspect  terrains  that  have  been  accreted  to 
the  eastern  North  American  continent  because  of  continental  collisions  during 
the  Paleozoic.  The  mechanism  of  plate  tectonics,  the  movements  and 
interactions  of  the  plates  that  form  the  earth’s  crust,  resulted  in  the 
opening  of  a  Late  Precambrian  to  Early  Paleozoic  lapetus  Ocean,  the  expansion 
of  the  lapetus  ocean  and  associated  deposition  of  sediments  into  this  ocean, 
and  the  eventual  closure  of  the  lapetus  ocean  by  the  Late  Paleozoic.  Times  of 
major  deformation  or  mountain  building  during  the  Paleozoic  are  interpreted  to 
correspond  to  periods  when  plate  collisions  have  occurred  and  crustal 
fragments  were  added  to  the  leading  edge  of  the  ancestral  North  American 
continent.  Three  main  periods  of  deformation  are  recognized  for  the 
Appalachian  Mountains.  These  periods  of  deformation  or  mountain  building 
correspond  approximately  to  the  Taconic  (450  to  500  m.y.),  Acadian  (350  to  400 
m.y.),  and  Alleghany  (250  to  300  m.y.)  Orogenies. 

10.  The  style  and  characteristics  of  deformation  during  each  tectonic 
period  varied  widely  along  the  Appalachian  chain  and  each  episode  of  mountain 
building  affected  segments  of  the  Appalachian  Mountains  differently.  Williams 
and  Hatcher  (1982)  interpret  the  Valley  and  Ridge,  the  Blue  Ridge,  and  the 
Piedmont  Provinces  as  being  suspect  terrains  that  approximately  correspond  to 
the  Piedmont,  Avalon,  and  Brunswick  terrains  as  shown  by  Figure  3  (from 
Wheeler  and  Bollinger,  1984).  Wheeler  and  Bollinger  (1984;  also  Bollinger  and 
Wheeler,  1983)  suggest  that  suspect  terrains  are  candidates  for  earthquake 
source  zones  and  that  they  may  help  explain  seismicity  in  the  southeastern 
United  States.  Causes  of  seismicity  for  the  southeastern  United  States  will 
be  examined  in  greater  detail  in  Part  III  of  this  report. 

11.  The  three  major  periods  of  deformation  during  the  Paleozoic 
produced  wide-spread  volcanism,  metamorphism,  folding,  and/or  large  scale 
westward  transport  of  numerous  vertically  stacked  thrust  sheets.  The  major 
period  of  thrust  faulting  and  folding  in  the  Central  and  Southern  Appalachians 
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Figure  2 .  Tectonic  map  of  Virginia  and  West  Virginia  showing  physiographic 
subdivisions  and  the  locations  of  major  faults  (fault  data  from  Price,  1986; 
Calver,  1963;  Butts,  1933;  Lowery  and  others,  1971;  and  Johnson,  1977) 


Fieure  3.  Map  showing  geologic  provinces  and  suspect  terrains  in  the 
southeastern  United  States  (from  Wheeler  and  Bollinger,  1984)  Provinces: 
small  dots  represents  Valley  and  Ridge;  stippled  represents  Blue  Ridge,  large 
dots  represents  Piedmont;  and  no  pattern  represents  the  Coastal  Plain.  Suspect 
terfain  boundaries  are  from  Williams  and  Hatcher  (1982).  Appa  ac  ian 
continental  margin  (AM),  Piedmont  (P) ,  Avalon  (A),  and  Brunswick  (B)  terrains 


is  interpreted  to  have  occurred  during  the  last  collision  event,  when  Africa 
and  North  America  were  joined  together  (Cook,  Brown,  and  Oliver,  1982;  Evans, 
1989;  Hatcher,  1972;  Lowery  and  others,  1971;  and  Van  der  Voo,  1979).  Thrust 
faulting  was  a  primary  mechanism  for  creation  of  the  Central  and  Southern 
Appalachian  Mountains  as  indicated  by  several  deep  penetrating  seismic 
reflection  profiles  (Harris,  deWitt,  and  Bayer,  1986;  Cook  Brown  and  Oliver, 
1981;  and  Oliver,  1982).  Underlying  the  Valley  and  Ridge  in  the  Central  and 
Southern  Appalachian  chain  is  the  interpreted  edge  of  the  proto-North  American 
continent.  Overlying  this  ancient  continental  edge,  are  the  thrust  faulted 
nearshore  and  marine  sediments  from  the  earlier  lapetus  Ocean.  These  thrusted 
sediments  form  the  underlying  geologic  section  at  Gathright  Dam.  West  of  the 
study  area  in  the  Appalachian  Plateau,  major  thrust  faults  are  absent; 
instead,  compression  was  limited  to  minor  folding. 

12.  The  final  stage  in  the  tectonic  history  of  the  Central  Appalachians 
and  the  southeastern  United  States  began  in  the  Mesozoic  Era  (250  to  65  m.y.). 
Rifting  separated  North  America  from  Africa  and  created  the  modern  Atlantic 
Ocean.  The  separation  of  the  two  land  masses  represents  a  change  in  the  style 
of  tectonism  from  compression  to  extension.  Relaxation  of  crustal  stresses 
produced  the  Triassic  basins  (250  to  210  m.y.),  which  are  bounded  by  normal 
faults,  and  produced  the  intrusion  of  numerous  cross-cutting  mafic  dikes 
throughout  the  southeast. 

13.  The  Triassic  basins  have  since  been  filled  with  sedimentary 
deposits  that  were  eroded  from  steep  mountains  to  the  east.  These  basins  are 
presently  buried  by  Cenozoic  deposits  (65  m.y.  to  present)  beneath  the  Coastal 
Plain.  Some  of  these  basins  are  exposed  in  the  Piedmont  Province  in  Virginia 
(Marine  and  Siple,  1974).  The  Triassic  mafic  dikes  are  well  developed  in  the 
southeastern  United  States  and  in  portions  of  central  and  northern  Virginia, 
where  they  cut  across  the  structural  grain  of  the  Appalachians,  extending 
approximately  northwest  to  southeast  (King  and  Beikman,  1976  and  Calver, 

1963).  Basin  formation,  normal  faulting,  and  dike  intrusion  are  interpreted 
to  have  ended  by  the  latter  part  of  the  Jurassic  Period  (210  to  145  m.y.  ago). 

14.  The  Cenozoic  Era  (65  m.y.  ago  to  present)  is  a  period  of  relative 
continental  stability.  The  coastal  plain  was  formed  during  this  time  from 
sediments  that  were  eroded  from  the  uplifted  Appalachian  Mountains  and 
deposited  along  the  continental  margin.  The  glacial  advances  during  the 
Pleistocene  (2  m.y,  to  10,000  years)  are  the  last  major  crustal  disturbances 
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to  have  occurred  in  North  America.  The  glaciers  did  not  advance  into  the 
Valley  and  Ridge  Province  (Flint,  1971). 


Regional  Geology  and  Structure 

15.  A  detailed  discussion  of  the  geology  in  the  Central  Appalachians  is 
beyond  the  scope  of  this  study.  The  following  discussion  will  be  restricted 
to  describing  the  fundamental  geologic  and  structural  characteristics  of  the 
Valley  and  Ridge  Province  as  they  relate  to  Gathright  Dam  and  to  possible 
sources  for  earthquakes. 

16.  The  major  geological  and  structural  features  in  Bath  and  Alleghany 
Counties  are  presented  in  Figure  4A  (from  Rader  and  Gathright,  1984)  with  the 
individual  stratigraphic  units  identified  in  Figure  4B  (from  Rader  and 
Gathright,  1984).  The  dam  and  reservoir  are  located  west  of  the  Warms  Springs 
Anticline.  The  geology  has  been  mapped  as  being  of  Paleozoic  age,  composed  of 
Ordovician  to  Mississippian  age  sediments.  Anticlines  and  synclines  are  the 
major  structural  features  in  the  three  county  area  identified  in  Figure  4A. 
This  area  of  Virginia  is  part  of  the  Western  Anticlines  (Rader  and  Gathright, 
1984).  Because  of  the  intense  folding  throughout  this  area,  the  structure  and 
stratigraphy  are  highly  variable.  The  sedimentary  rocks  that  comprise  the 
area  have  been  subjected  to  several  periods  of  deformation,  producing  multiple 
fold  orientations  with  plunging  anticlinal  and  synclinal  structures. 

17.  Bath  and  Highland  Counties  are  both  noted  for  unusual  geologic 
features.  These  two  counties  contain  the  largest  concentration  of  thermal 
springs  in  the  eastern  United  States  (Bollinger  and  Gilbert,  1974).  These 
thermal  springs  are  often  associated  with  travertine  deposits.  The  locations 
of  hot  springs  in  close  proximity  to  Lake  Moomaw  and  Gathright  Dam  are 
identified  in  Figure  1  (note  the  town  locations).  In  addition.  Highland 
County  contains  Eocene  age  (58  to  37  m.y.)  volcanic  intrusions,  the  youngest 
known  igneous  rocks  in  the  eastern  United  States  (Dennison  and  Johnson,  1971; 
and  Kettren,  1971).  Dennison  and  Johnson  (1971)  have  proposed  that  a  cooling 
igneous  body  in  the  subsurface  is  responsible  for  the  hot  springs  and  once  was 
the  source  for  the  igneous  intrusions  in  Highland  County.  The  existence  of 
the  igneous  intrusion  has  yet  to  be  proven.  Studies  by  Costain  and  others 
(1976  and  1978)  do  not  support  the  existence  of  the  cooling  intrusion; 
instead,  they  suggest  an  alternative  hypothesis  of  deep  circulation  of  surface 
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Figure  4B.  Stratigraphic  units  in  the  Western  Anticlines  region  and  at 
Gathright  Damsite  (from  Rader  and  Gathright,  1984) 
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water  on  a  normal  geothermal  gradient  by  way  of  existing  fractures  and  faults. 

18.  The  major  structures  in  the  Western  Anticlines  area  are  related  to 
the  compressional  events  that  occurred  during  the  Paleozoic  when  major  thrust 
faulting  and  folding  are  interpreted  to  have  occurred.  A  generalized 
subsurface  section  of  the  Western  Anticlines  area,  approximately  along  the 
Alleghany  and  Bath  County  line  (see  Figure  4A) ,  is  presented  in  Figure  5 
(modified  from  Kulander  and  Dean,  1978).^  The  structural  interpretations 
shown  in  Figure  5  are  based  on  detailed  gravity  and  magnetic  surveys  along  the 
Alleghany  Plateau  and  portions  of  the  Valley  and  Ridge  Province  near  the  Warm 
Springs  Anticline.  The  Warm  Springs  Anticline  is  interpreted  to  be  a  major 
stratigraphic  flexure  that  is  part  of  and  related  to  a  series  of  deep  seated 
thrust  faults  that  sole  at  depth  into  a  master  decollement  surface  above  the 
Precambrian  surface. 

19.  The  local  stratigraphy  and  structure  along  the  Valley  and  Ridge  and 
Appalachian  Plateau  boundary  will  vary  with  location  along  the  boundary,  but 
the  overall  mechanism  of  compressional  folding  and  thrust  faulting  is  common 
along  the  entire  length  for  both  the  southern  and  central  segments  of  the 
Appalachian  Mountain  chain.  However;  the  major  linear  faults  so 
characteristic  of  the  Southern  Appalachians  (see  Figure  2)  are  not  continuous 
into  the  Central  Appalachians.  The  presence  of  thrust  faults  at  the  surface 
diminishes  in  the  central  segment  as  compared  to  the  southern  segment.  The 
presence  of  major  thrust  faults  in  the  subsurface  is  well  noted  for  the 
central  segment  (Evans,  1989;  Wilson,  1989;  Kulander  and  Dean,  1978;  and  Dean 
and  Kulander,  1972).  Further  to  the  west  in  the  Appalachian  Plateau  region, 
the  effects  of  major  thrust  faulting,  common  in  the  Valley  and  Ridge,  the  Blue 
Ridge,  and  the  Piedmont  Provinces,  are  negligible.  Instead,  Paleozoic 
compression  in  the  sedimentary  cover  is  expressed  by  minor  folding  that 
eventually  diminishes  and  merges  with  the  relatively  flat- lying  depositional 
sequences  of  the  central  continent.  The  deformation  of  the  sedimentary  cover 
in  the  central  Appalachians  and  in  the  project  area  is  interpreted  as  being 
one  of  thin  skinned  tectonism.  deformation  of  the  upper  crust  without  active 
involvement  of  the  underlying  crystalline  basement  rocks  (Rodgers,  1971  and 


^The  cross-section  in  Figure  5  is  modified  from  the  original  northwest 
to  southeast  cross-section  K  -  K'  by  Butts  (1933),  located  approximately 
1.5  miles  south  of  the  damsite. 
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Figure  5.  Generalized  subsurface  structure  of  Western  Valley  and  Ridge 
Province  near  vicinity  of  Alleghany  and  Bath  County  line,  Virginia.  Subsur  ace 
structure  interpreted  from  detailed  gravity  and  magnetic  data  (from  Kulander 
and  Dean,  1978;  cross-section  location  and  geology  modified  from  cross-section 

K-K'  by  Butts,  1933) 
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Gwinn,  1964). 


General  Site  Geology  and  Structure 

20.  General  information  concerning  the  site  geology  and  the  structural 
features  in  the  project  area  are  summarized  below  and  were  obtained  from 
various  design  memoranda  and  foundation  reports  (U.S.  Army  Corps  of  Engineers, 
1966,  1967,  1969,  1976,  1983a,  and  1983b).  The  primary  geologic  units  and 
structural  features  underlying  Gathright  Dam  and  Lake  Moomaw  are  identified  in 
Figure  6  (from  U.S.  Army  Corps  of  Engineers,  1983b). 

21.  Gathright  Dam  and  Lake  Moomaw  are  underlain  by  Silurian  (438  to 
408  m.y.)  and  Devonian  (408  to  350  m.y.)  sedimentary  rocks  that  have  been 
folded  and  faulted.  The  major  structural  features  in  the  dam  and  reservoir 
area  are  folds,  faults,  and  joints.  As  previously  noted,  the  intense  folding 
during  the  Paleozoic  produced  anticlinal  and  synclinal  structures.  Bolar 
Mountain,  Coles  Mountain,  Morris  Hill,  and  Hoover  Ridge  are  anticlines. 
Gathright  Dam  is  located  in  a  narrow  gorge  which  the  Jackson  River  has  cut 
between  Bolar  and  Coles  Mountains.  The  dam  is  located  a  short  distance 
upstream  from  the  axis  of  the  Morris  Hill  anticline.  Lake  Moomaw  is  contained 
within  a  syncline  between  the  surrounding  anticlinal  structures.  Detailed 
information  on  the  site  geology  and  important  structural  features  are 
presented  in  Appendix  A.  Included  in  Appendix  A  is  a  geologic  cross-section 
from  beneath  the  dam  that  shows  the  orientation  and  distribution  of  the 
various  rock  units. 


Determination  of  Active  and  Capable  Faults 

Definition  of  Active  and  Capable  Fault 

22.  Earthquakes  are  produced  when  strain  energy  is  suddenly  released  in 
the  form  of  movements  along  faults.  The  identification  and  recognition  of 
active  faults  are  important  in  determining  the  earthquake  potential  for  an 
area.  An  active  fault  is  defined  by  the  U.S.  Army  Corps  of  Engineers  (1983c) 
as  a  fault  which  has  moved  during  the  recent  geologic  past  (Quaternary)  and, 
consequently  may  move  again.  However,  an  active  fault  may  or  may  not  be 
capable  of  producing  earthquakes.  An  active  fault  is  judged  capable  of 
producing  earthquakes  if  it  is  shown  to  exhibit  one  or  more  of  the  following 
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Figure  6.  Geology  and  structure  at  Gathright  Dam  and  vicinity  (from  U.S.  Army 
Corps  of  Engineers,  1983b;  Plate  II-3).  Detailed  descriptions  of  the 
individual  rock  units  are  presented  in  Appendix  A 
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characteri53tics ; 


a.  Movement  at  or  near  the  ground  surface  at  least  once  during  the  past 
35,000  years. 

b.  Macroseismicity  (M^3.5)  instrumentally  determined  with  records  of 
sufficient  precision  to  demonstrate  a  direct  relationship  with  the 
fault . 

c.  A  structural  relationship  to  a  capable  fault  such  that  movement  on 
one  fault  could  be  reasonably  expected  to  cause  movement  on  the  other. 

d.  Established  patterns  of  microseismicity  that  define  a  fault  and 
historic  macroseismicity  that  can  reasonably  be  associated  with  that 
fault . 

23.  In  summary,  a  fault  that  is  active  and  capable  of  producing 
earthquakes  must  show  either  geologic  and/or  seismic  evidence  of  its  activity. 
Geologic  Evidence 

24.  Numerous  geologic  studies  have  been  conducted  in  the  central 
Appalachian  region  to  identify  lineaments  and  possible  Post-Paleozoic  fault 
activity  (White,  1952;  Mixon  and  Newell,  1977;  Dames  and  Moore,  1977; 

Reynolds,  1979;  Geiser,  1978;  Wheeler,  1980;  Rader  and  Gathright,  1984;  and 
Southworth,  1986a  and  1986b).  These  studies  are  mostly  of  a  general  or 
reconnaissance  nature  that  describe  major  lithologic  aspects  of  the  different 
tectonic  events,  significant  structural  characteristics  of  ancient  tectonism, 
interpretations  of  basement  structure,  or  structural  characteristics  that  are 
related  to  mineral  and  petroleum  exploration.  For  the  most  part,  these 
studies  do  not  identify  known  Cenozoic  age  faults  or  possible  earthquake 
sources  at  or  near  Gathright  Dam. 

25.  One  study  that  reports  a  possible  Tertiary  age  fault  near  Gathright 
Dam  is  by  White  (1952).  Faulting  was  identified  in  a  road  cut  on  U.S.  Highway 
60  near  Clifton  Forge.  The  road  cut  is  approximately  1  mile  due  east  of  the 
junction  between  U.S.  Highways  220  and  60,  approximately  12  miles  (19  km) 
southeast  from  Gathright  Dam.  The  fault  is  described  by  White  (1952)  as  being 
vertical.  It  displaces  a  surface  gravel  horizon  that  was  fluvially  deposited 
and  an  underlying  shale  bed.  The  fault  strike  is  N  25°  W,  or  transverse  to 
the  strike  of  the  central  Appalachian  Mountain  chain.  White  indicates  that 
the  fault  is  of  tectonic  origin,  based  on  geomorphic  evidence,  and  he  assumed 
it  to  be  Tertiary  in  age.  Other  examples  of  normal  and  reverse  Cenozoic 


faulting  are  cited  by  White  for  several  other  locations  in  Virginia  and  North 
Carolina.  He  concludes  that  many  more  such  minor  faults  may  be  present  in  the 
southeast.  To  date,  there  are  no  proven  Holocene  (less  than  10,000  years) 
faults  identified  in  the  literature  for  the  southeastern  United  States. 

26.  The  U.S.  Army  Corps  of  Engineers,  Norfolk  District,  conducted  a 
detailed  evaluation  of  faulting  as  part  of  the  design  and  construction  of 
Gathright  Dam.  As  part  of  the  evaluation  process,  the  Norfolk  District 
conducted  an  extensive  review  of  the  literature  and  performed  an  in-depth 
geological  analysis  of  the  impoundment  area  in  1973.  The  geological  analysis 
was  conducted  by  the  Raytheon  Company  under  the  direction  of  the  Office,  Chief 
of  Engineers.  The  analysis  by  Raytheon  used  color  aerial  photography  and  high 
resolution  classified  black  and  white  imagery. 

27.  The  faults  and  lineaments  that  were  identified  by  the 
Norfolk/Raytheon  study  were  mapped  on  approximately  1:24,000  scale  base  maps 
and  these  features  were  then  evaluated  in  the  field  by  a  geologist.  The 
mapped  faults  in  the  Gathright  Dam  impoundment  area  and  vicinity  are  shown  in 
Figure  6.  The  final  report  that  describes  the  work  performed  and  results 
obtained  from  this  study  on  faulting  was  classified.  There  are  presently  no 
copies  of  this  report  (U.S.  Army  Corps  of  Engineers,  1973)  available  as  it  was 
destroyed  for  security  reasons.  However,  results  of  this  detailed  study  have 
been  summarized  in  two  later  Corps  reports  (U.S.  Army  Corps  of  Engineers,  1976 
and  1983b).  It  was  concluded  by  the  Norfolk  District,  based  on  their  in-depth 
studies,  that  there  were  no  active  and  capable  faults  in  the  project  area. 

28.  The  following  activities  were  performed  as  part  of  this  study  to 
identify  and  evaluate  faults  near  Gathright  Dam  and  to  determine  whether  any 
of  these  faults  are  active  and  capable  of  producing  earthquakes. 

a.  An  extensive  review  of  the  literature  was  conducted  to  evaluate  the 
more  recent  geologic  studies  for  evidence  of  tectonism,  seismicity,  and 
the  presence  of  active  faults  in  the  southeastern  United  States. 

b.  Aerial  photography  (black  and  white  1:24,000,  black  and  white 
1:20,000,  and  color  infrared  1,24,000  scale)  were  examined  to  identify 
faults  and  linears  at  Gathright  Dam  and  the  surrounding  vicinity. 

c.  Obtain  and  review  technical  information  and  interpretations  from 
government  and  university  geologists  and  seismologists  who  have 
knowledge  about  the  geology  and  seismicity  of  the  study  area. 
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d.  Review  and  evaluate  the  historic  record  of  seismicity  in  order  to 
establish  causative  relations  with  known  faults.  A  review  and  analysis 
of  the  seismic  history  and  seismicity  for  the  study  area  is  presented  in 
the  next  section  of  this  report. 

29.  There  is  no  evidence  in  the  surface  geology  or  the  seismic  history 
that  identifies  faults  capable  of  producing  earthquakes  at  or  near  Gathright 
Dam.  It  is  concluded  after  evaluating  the  above  information  that  there  are  no 
identifiable  active  or  capable  faults  at  or  near  Gathright  Dam. 
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PART  III:  SEISMICITY 


Relation  Between  Seismicity  and  Geology 

30.  Geophysical  studies  are  useful  in  identifying  anomalous  geologic 
structures  deep  within  the  subsurface.  These  structures  may  indicate  areas 
where  tectonic  stresses  are  concentrated  and  whose  potential  sources  exist  for 
earthquakes.  The  sudden  release  of  built-up  strain  energy  from  the  focusing 
of  tectonic  stresses  produces  earthquakes.  Gravity  and  magnetic  surveys  are 
two  important  types  of  geophysical  studies  that  help  to  define  geological 
irregularities  or  structures  in  the  subsurface  where  regional  tectonic 
stresses  may  concentrate. 

31.  Figure  7  presents  the  results  of  a  generalized  gravity  survey  for 
Virginia  and  West  Virginia  (modified  from  Johnson,  1977;  and  Dean  and 
Kulander,  1987).  A  gravity  map  identifies  density  variations  which  in  turn 
indicate  differences  in  rock  type  and  thickness.  The  gravity  map  closely 
parallels  the  trend  of  the  Appalachian  Mountains.  It  defines  the  general 
boundaries  between  the  different  tectonic  provinces  by  variations  in  the 
contour  intensity.  The  boundary  between  the  Valley  and  Ridge  and  Blue  Ridge 
provinces  is  defined  by  the  steep  gravity  gradient  beginning  at  the  Blue  Ridge 
and  Valley  and  Ridge  Boundary.  To  the  east  of  the  Valley  and  Ridge  province, 
the  gravity  map  defines  a  gravity  high  that  extends  to  the  coastal  plain.  To 
the  west  of  the  Valley  and  Ridge  province,  the  gravity  map  detines  a  broad 
gravity  low.  Gathright  Dam  is  centered  over  the  deepest  part  of  the  gravity 
low. 

32.  Kulander  and  Dean  (1978)  identify  the  gravity  low  beneath  Gathright 
Dam  as  corresponding  to  the  thickest  section  of  crust  in  the  central 
Appalachians  with  an  estimated  thickness  of  approximately  34  miles  (55  km). 
They  indicate  that  the  low  beneath  the  Gathright  Dam  represents  the 
combination  of  the  sum  gravational  effect  of:  a)  large  crustal  thickness  and 
lateral  compositional  variations  in  crust  and  mantle,  b)  broad  basement 
surface  with  low  relief,  and  c)  tectonically  produced  density  variations  in 
the  cover  of  Paleozoic  sedimentary  rocks. 

33.  An  aeroraagnetic  map  is  presented  in  Figure  8.  The  aeromagnetic  map 
identifies  areas  having  a  susceptibility  or  remnant  magnetization  of 
sufficient  magnitude  to  produce  a  measurable  distortion  in  the  earth's 
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magnetic  field.  Igneous  rocks  are  the  primary  sources  for  magnetic  minerals 
that  are  capable  of  producing  variations  in  the  magnetic  field.  The 
aeromagnetic  map  in  general  identifies  the  structural  orientation  of  the 
central  Appalachians.  Highest  values  occur  in  the  Blue  Ridge  province  where 
the  bulk  of  the  igneous  rocks  are  concentrated.  Beneath  Gathright  dam  there 
are  no  anomalous  magnetic  structures.  The  absence  of  a  sharp  magnetic  anomaly 
is  due  to  the  thick  sequence  of  thrust  faulted  sedimentary  rocks  which  occur 
at  this  location  (see  Figure  5) .  The  aeromagnetic  map  in  general  corroborates 
the  boundaries  and  other  tectonic  features  identified  by  the  gravity  map  and 
indicates  an  absence  of  anomalous  structures  near  the  damsite. 

Causes  of  Earthquakes 

34.  Earthquakes  are  produced  when  strain  energy  is  suddenly  released  in 
the  form  of  movements  along  faults.  Strain  energy  is  derived  from  regional 
tectonic  stresses  which  are  created  by  the  interactions  between  the  major 
crustal  plates  that  form  the  earth's  surface.  The  sudden  movement  or  slippage 
along  fault  surfaces  produces  an  elastic  rebound.  This  elastic  rebound 
generates  vibrations  in  the  earth's  crust  and  these  vibrations  are  felt  as  an 
earthquake.  To  have  a  large  earthquake  requires  a  large,  sudden  stress  drop 
and  energy  release  which  can  only  be  produced  by  fault  movements  that  are 
rooted  in  the  deep  crystalline  basement  rocks. 

35.  The  causes  of  earthquakes  in  the  southeastern  United  States  and  the 
study  area  are  not  well  understood  as  active  and  capable  faults  have  not  been 
identified  to  date  at  the  surface.  Since  active  faults  have  not  been 
identified  at  surface,  there  are  several  theories  that  have  been  proposed  to 
explain  the  causes  of  earthquakes  in  the  southeastern  United  States: 

a.  Focusing  of  regional  stresses  at  heterogeneities  (plutons)  or  other 
discordant  rock  masses  in  the  subsurface  and  release  of  this  stress  by 
fault  movements  at  depth. 

b.  Introduction  of  small-scale  magmatic  materials  into  the  lower  crust, 
producing  stresses,  and  generating  fault  movements  at  depth. 

c.  Focusing  and  release  of  regional  stresses  along  major  tectonic 
discontinuities  such  as  ancient  rift  zones  or  transform  faults. 

d.  Regional  compression  causing  activation  and  slippage  along 
preexisting  faults  planes  such  as  thrust  faults. 
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e.  Regional  extension  producing  movements  along  fault  bounded  coastal 
graben  structures  (Triassic  Basins)  or  relaxation  type  movements  on 
existing  faults  (Barosh,  1981  and  1987;  and  Armbruster  and  Seeber, 

1981) . 

f.  Localized  stress  relief  along  joint  planes  or  other  near  surface 
discontinuities  (Long,  1988;  and  Costain,  Bollinger,  and  Speer,  1987). 
These  earthquakes  are  related  to  ground  water  movements  and  water  table 
fluctuations . 

36.  The  generally  accepted  view  of  eastern  United  States  seismicity  is 
that  earthquakes  occur  along  pre-existing  zones  of  weakness  that  are  favorably 
orientated  with  respect  to  the  present  stress  field  (Bollinger,  Ehlers ,  and 
Moses,  1987).  These  pre-existing  zones  of  weakness  are  interpreted  to  be 
ancient  faults,  paleorift  zones,  or  the  intersection  points  for  multiple 
tectonic  features  which  are  located  deep  in  the  subsurface.  Bollinger, 

Ehlers,  and  Moses  (1987)  indicate  the  major  concentration  of  earthquakes  in 
the  southeastern  United  States  occurs  along  the  Avalon  and  Piedmont  boundary 
(see  Figure  3).  Barosh  (1987)  points  out  that  all  of  the  major  seismic  areas 
in  the  eastern  United  States  lie  in  or  adjacent  to  the  heads  of  Late 
Cretaceous-Early  Tertiary  embayments.  In  areas  where  these  Cretaceous-Early 
Tertiary  embayment  type  deposits  are  absent  and  historic  earthquakes  have 
occurred,  Barosh  attributes  movements  along  pre-existing  zones  of  weakness, 
either  to  sedimentary  loading  related  to  these  basins  or  to  fault  reactivation 
related  to  the  formation  of  the  Atlantic  Ocean  and  the  Gulf  of  Mexico. 

37.  The  exact  cause  of  seismicity  in  the  eastern  United  States  will  be 
dependent  on  the  site  geology  and  defined  by  the  historic  record  of 
seismicity.  Explanations  a  through  e  in  the  above  list  can  be  interpreted  as 
suggesting  that  large  earthquakes  can  happen  anywhere  in  the  eastern  United 
States  or  the  study  area  at  a  location  where  no  historic  earthquakes  have 
occurred  in  the  past.  To  project  an  earthquake  into  an  area  or  a  zone  that 
has  displayed  no  past  seismicity,  but  is  part  of  a  major  ancient  fault  trend 
or  zone,  is  not  considered  valid  by  the  present  authors  unless  there  is 
additional  evidence  for  seismicity.  A  key  question  that  must  be  asked  in  such 
an  evaluation  as  this,  is  there  a  relation  between  the  present  seismicity  and 
the  existing  geologic  structures?  The  folding  and  faulting  that  have  been 
mapped  for  the  central  Appalachians  (see  Figures  2  and  6)  are  derived  from 
ancient  tectonism  which  is  no  longer  active  today.  Present  day  tectonism  is 
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greatly  different  from  the  tectonisra  which  formed  these  ancient  structures. 

The  present  seismicity  is  related  to  the  tectonism  and  stress  fields  which  are 
active  today. 

38.  Explanation  f  above  implies  a  very  low  upper  bound  on  the  maximum 
earthquake  that  can  occur.  The  maximum  intensity  level  for  this  upper  bound 
is  believed  to  be  at  a  level  that  is  below  that  of  concern  for  engineering. 
Stress  release  is  near  the  surface,  generally  unrelated  to  any  geologic 
structures  except  for  joints.  Some  earthquakes  are  believed  to  be  triggered 
by  ground-water  movements  through  the  joints.  However,  since  these 
earthquakes  are  shallow  and  of  low  energy,  a  major  earthquake  is  not  expected 
to  be  generated  by  this  mechanism.  In  addition,  this  type  of  earthquake  would 
be  especially  apt  to  occur  near  reservoirs. 

39.  Microearthquake  monitoring  was  conducted  as  part  of  the 
construction  of  the  Bath  County  Pump  and  Storage  facility,  a  major 
hydroelectric  facility  in  northern  Bath  County  (approximately  10  miles,  16  km) 
north  of  Gathright  Dam.  The  Bath  County  Pump  and  Storage  project  included  the 
construction  of  two  large  storage  reservoirs  (see  Figures  1  and  9  for  location 
of  Back  Creek  and  Little  Back  Creek  Reservoirs) .  Microearthquake  monitoring 
was  discontinued  after  several  years  following  construction  of  the 
hydroelectric  facility  as  there  was  no  significant  earthquake  activity 
associated  with  the  filling  of  the  reservoirs  or  the  daily  discharge  and 
filling  of  the  two  reservoirs  (Bollinger,  personnel  communication). 

Information  on  microearthquake  monitoring  in  the  project  area  and  at  the  Bath 
County  hydroelectric  facility  is  presented  in  a  later  section  of  this  report. 
It  is  concluded  that  earthquakes  from  "hydroseismic"  sources  will  not  produce 
large  earthquakes  in  this  region  that  would  adversely  affect  engineered 
structures . 

40.  In  summary,  the  maximum  earthquake  potential  is  a  function  of  the 
present  day  tectonism.  It  must  be  assumed  that  the  largest  earthquakes  that 
can  occur  in  the  study  area  are  defined  by  the  seismic  history  and/or  by  the 
presence  of  active  faults.  These  two  considerations  will  control  the 
selection  process  for  the  maximum  earthquake  that  is  specified. 
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Al.  The  distribution  of  historic  earthquakes  in  the  study  area  is  pre¬ 
sented  in  Figure  9.  Earthquakes  are  shown  according  to  the  Modified  Mercalli 
Intensity  (MMI)  scale.  The  catalogue  of  historic  earthquakes  from  which 
Figure  9  was  prepared  is  contained  in  Appendix  B.  The  catalogue  is  derived 
from  the  Earthquake  Data  Base  maintained  by  the  National  Geophysical  Data 
Center,  National  Atmospheric  and  Oceanic  Administration,  Boulder,  Colorado 
(from  Habermann,  1989) .  The  list  of  historic  earthquakes  is  arranged  by  date 
and  time  (Universal  or  Greenwich  Time)  and  includes  coordinate  location  of  the 
epicenter,  earthquake  magnitude  (mj,.  Ml,  and  M^) ,  MMI,  and  focal  depth.  A 
glossary  of  terms  is  included  in  Appendix  C  which  includes  a  description  of 
the  MMI  scale  and  the  different  instrumental  or  magnitude  scales  that  are 
used. 

42.  The  catalogue  in  Appendix  B  contains  a  listing  of  46  events  between 
the  years  18ol  and  1986,  a  history  that  spans  185  years.  The  vast  majority  of 
earthquakes  are  less  than  MMI  IV.  Only  four  earthquakes  are  larger  than  MMI 
V.  The  largest  earthquake  in  the  catalogue  is  MMI  VIII.  The  distribution  of 
historic  earthquakes  is  as  follows;  2  earthquakes  of  MMI  I,  5  earthquakes  of 
MMI  II,  8  earthquakes  of  MMI  III,  10  earthquakes  of  MMI  IV,  14  earthquake  of 
MMI  V,  3  earthquakes  of  MMI  VI,  and  1  earthquake  of  MMI  VIII.  The  MMI  VIII 
earthquake  occurred  on  31  May  1897  in  Giles  County,  Virginia. 

43.  Examination  of  Figure  9  and  the  State  Seismicity  Map  of  Virginia 
(Reagor,  Stover,  and  Algermissen,  1983)  indicates  that  there  is  no  significant 
concentration  of  historic  earthquakes  in  the  study  area  other  than  at  Giles 
County.  The  highest  concentration  of  earthquakes  in  the  project  area  occurs 
at  Giles  County.  There  are  no  significant  historic  earthquakes  located  at  or 
near  Gathright  Dam.  The  seismic  record  indicates  a  region  surrounding 
Gathright  Dam  which  is  classified  as  aseismic  based  on  the  distribution  of 
historic  earthquakes. 


Microear thouakes 


Bath  County  Microearthauakes 

44.  Microearthquakes  are  earthquakes  that  are  too  small  to  be  felt 
(i.e.,  M  ^3),  but  are  recorded  by  seismographic  instruments. 
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Microearthquakes  are  useful  for  defining  areas  where  tec  nic  stresses  are 
concentrated.  These  small  earthquakes  are  helpful  in  determining  focal 
depths,  fault  types  and  their  orientations,  and  aid  in  estimating  rates  of 
earthquake  recurrence.  Microearthquakes  are  important  in  determining  if  there 
is  a  correlation  between  ancient  tectonic  structures  (i.e.,  faults,  plutons, 
etc.)  and  present  day  seismic  activity. 

45.  Microearthquake  monitoring  has  been  conducted  in  Bath  County  as 
part  of  the  construction  and  seismic  evaluation  of  Virginia  Power  and  Electric 
Company's  new  Bath  County  Pump  and  Storage  Facility,  a  2.1  megawatt 
hydroelectric  plant.  The  Bath  County  facility,  located  approximately  10  miles 
(16  km)  north  of  Gathright  Dam,  consists  of  two  reservoirs  (see  Figure  1  for 
location  of  Little  Back  Creek  and  Back  Creek  Reservoirs),  a  powerhouse,  and 
related  facilities  necessary  for  generating  hydroelectric  power.  A  crustal 
velocity  model  was  developed  for  Bath  County  from  construction  blasts  to 
accurately  define  earthquake  locations  and  focal  depths.  Microearthquake 
monitoring  was  conducted  for  a  3  month  period  during  1982  with  a  portable 
seismographic  network  and  also  beginning  in  1978  with  a  permanent  4  station 
network  (Todd,  1982).  Todd  (1982)  concluded  that  Bath  County  was  aseismic. 

The  three  month  monitoring  program  failed  to  detect  any  local  earthquake 
activity.  Only  eleven  earthquakes  were  reported  for  the  long-term  monitoring 
program  between  1978  and  1982  and  three  of  these  events  were  to  small  to 
locate.  The  permanent  monitoring  program  was  finally  discontinued  in  1987  as 
there  was  very  little  microearthquake  activity.  There  was  no  activity 
associated  with  the  initial  filling  of  the  two  reservoirs  or  from  the  daily 
filling  and  emptying  of  these  reservoirs  as  part  of  the  hydroelectric 
generating  process  (Bollinger,  personnel  communication). 

46.  The  U.S.  Army  Corps  of  Engineers,  Norfolk  District,  did  not  conduct 
a  program  to  monitor  microearthquakes  during  filling  of  Lake  Moomaw.  However, 
the  Bath  County  Pump  and  Storage  Facility  and  the  Giles  County  seismic 
networks  were  operational  during  this  period  and  would  have  detected  any 
anomalous  activity  associated  with  the  reservoir  filling  at  Gathright  Dam.  To 
date,  there  have  been  no  published  reports  of  any  anomalous  microseismic 
activity  in  Bath  County,  associated  with  either  Gathright  Dam  or  the  Bath 
County  Pump  and  Storage  Facility. 

Virginia  Microearthauakes 

47.  Bollinger  and  others  (1986)  identify  two  source  areas  for 
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pronounced  microearthquake  activity  in  Virginia  as  shown  in  Figure  10.  The 
catalogue  from  which  Figure  10  is  derived  is  presented  in  Appendix  D  (from 
Bollinger  and  others,  1986).  The  major  source  of  microearthquake  activity  in 
the  project  area  occurs  in  Giles  County,  the  source  area  for  Virginia's 
largest  historic  earthquake  and  the  second  largest  historic  earthquake  in  the 
southeastern  United  States.  A  second  zone  of  microearthquake  activity  occurs 
outside  of  the  project  area,  within  Virginia's  Piedmont  Province,  and  is  known 
as  the  Central  Virginia  seismic  zone.  Earthquake  monitoring  in  Virginia 
between  1977  and  1985  indicates  that  significant  differences  occur  between  the 
two  seismic  source  areas  (Bollinger  and  others,  1986;  Bollinger  and  Wheeler, 
1982,  1982a,  and  1983;  and  Munsey,  1984).  In  Giles  County,  seismic  energy  is 
predominantly  released  by  strike  slip  movements  along  a  near  vertical,  tabular 
zone,  approximately  25  miles  (40  km)  long.  Movements  are  in  the  crystalline 
basement  rock,  approximately  3  to  16  miles  (5  to  25  km)  deep,  and  below  the 
basal  detachment  of  the  thrust  faulted  Appalachian  sediments  (see  Appendix  D 
for  focal  depths).  In  central  Virginia,  the  release  of  seismic  energy  is 
above  the  crystalline  basement  rock,  in  the  thick  stack  of  thrust  faulted 
Appalachian  sediments.  Movements  are  from  a  combination  of  dip  slip  and 
strike  slip  movements  that  occur  within  a  circular  area  approximately  62  miles 
(100  km)  in  diameter  and  6  miles  (10  km)  in  vertical  thickness. 

48.  A  comparison  of  microeatthquake  epicentral  locations  with  the 

epicentral  locations  for  historic  earthquakes  between  the  years  1774  to  1977, 
prior  to  instrumentation,  is  shown  by  Figure  11  (from  Bollinger  and  others, 
1986) .  The  source  areas  for  microearthquakes  are  in  close  correspondence  with 
the  locations  of  historic  earthquakes.  The  above  distribution  suggests  that 
earthquake  activity  in  Virginia  is  confined  to  "hot  spots"  which  are  defined 
by  present  day  microearthquake  activity.  Gathright  Dam  and  Lake  Moomaw  are 
situated  in  an  area  with  very  little  of  either  category  of  earthquakes.  There 
has  been  no  significant  concentration  of  historic  earthquakes  or  anomalous 
microearthquake  activity  at  the  darasite.  Gathright  Dam  appears  to  be  located 
in  an  aseismic  area.  The  major  source  for  earthquake  activity  in  the  project 
area  is  located  at  Giles  County.  * 

49.  Bollinger  and  Wheeler  (1982)  have  published  a  very  comprehensive 
evaluation  and  analysis  of  the  geology,  the  historic  seismicity,  the 
microcarthquake  activity,  and  the  tectonism  in  Giles  County.  Their 
interpretation  for  the  causes  of  earthquake  activity  at  Giles  County  is  that 
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Figure  10.  Distribution  of  instrumentally  located  earthquakes  in  Virginia. 
Earthquakes  shovm  in  this  figure  are  presented  in  Appendix  C  (from  Bollinger 

and  others,  1987) 


earthquakes  are  produced  from  compressional  reactivation  of  laoetan  normal 
faults  which  occur  in  the  basement  in  response  to  the  present  stress  field. 

Seismic  Source  Zones  in  the  Southeastern  United  States 

50.  Earthquake  source  zones  are  interpreted  for  the  southeastern  United 
States  as  there  are  no  known  active  faults.  These  source  zones  are  based  on 
the  record  of  historic  and  instrumentally  recorded  earthquakes.  The  seismic 
source  zones  interpreted  for  the  southeastern  United  States  are  shown  in 
Figure  12.  The  southeastern  United  States  is  in  general  a  region  of  low  level 
seismicity  with  areas  of  concentrated  earthquake  activity.  These  concentrated 
areas  or  zones  are  called  seismic  "hot  soots"  and  are  potential  sources  for 
moderate  to  major  earthquakes. 

51.  An  earthquake  zone,  as  used  in  this  report,  is  an  inclusive  area 
over  which  a  given  maximum  credible  earthquake  can  occur.  This  is  the  largest 
earthquake  that  can  reasonably  be  expected  to  occur.  It  can  be  moved  anywhere 
in  the  zone  and  is  thus  a  floating  earthquake.  The  earthquake  is  moved  in 
this  manner  because  causative  faults  have  not  been  identified.  The  criteria 
by  which  the  seismic  zones  in  Figure  12  were  developed  are  as  follows: 

a.  Maximum  sizes  of  earthquakes. 

b.  Density  of  earthquakes,  using  historic  seismicity  plus  microseismic 
activity  where  available.  A  strong  occurrence  of  both  together 
identifies  a  seismic  hotspot. 

c.  One  earthquake  will  adjust  a  boundary  but  cannot  create  a  zone. 

d.  Zones  of  greatest  activity  are  generally  as  small  as  possible. 

e.  The  maximum  intensity  of  a  zone  cannot  be  smaller  but  may  be 
equal  to  or  greater  than  the  maximum  historic  earthquake. 

f.  The  zones  are  source  areas.  They  do  not  necessarily  represent  the 
maximum  intensity  at  every  point  since  attenuations  have  to  be  taken 
into  account. 

52.  The  largest  earthquake  source  zones  in  this  portion  of  the  United 
States  are  at  Charleston,  South  Carolina  and  at  Giles  County,  The  Charleston 
area  is  shown  as  generating  an  earthquake  of  MMI  X.  An  intensity  MMI  X 
earthquake  occurred  at  Charleston  in  1886.  The  Giles  County  area  is  shown  as 
having  the  potential  to  generate  an  earthquake  of  MMI  IX.  The  Giles  County 
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Figure  12.  Seismic  source  zones  in  the  southeastern  United  States 


source  area  is  interpreted  as  being  one  intensity  level  higher  than  the 
largest  historic  earthquake  that  has  occurred.  As  previously  stated,  an 
earthquake  of  MMI  VIII  occurred  at  Giles  County  in  1897. 

53.  Examination  of  Figure  12  shows  Gathright  Dam  as  being  located  in  an 
area  experiencing  low  level  seismic  activity.  The  maximum  earthquake 
interpreted  for  the  damsite  location  is  MMI  VI.  Gathright  Dam  is  bordered  on 
the  southwest  by  the  Giles  County  seismic  zone  (MMI  IX) ,  on  the  east  by  the 
Central  Virginia  seismic  zone  (MMI  VII) ,  and  on  the  northeast  by  Northern 
Virginia  seismic  zone  (MMI  VII). 

54.  The  Giles  County  and  the  Central  Virginia  seismic  zones  are  both 
located  within  50  miles  (80  km)  of  Gathright  Dam.  The  Northern  Virginia 
Seismic  zone  is  located  more  than  50  miles  (80  km)  from  Gathright  Dam. 

Because  of  its  close  proximity  to  Gathright  Dam,  Giles  County  is  the 
controlling  earthquake  source  area  for  Gathright  Dam.  Consequently,  the 
Central  Virginia  and  Northern  Virginia  seismic  zones  are  not  considered  to  be 
as  important  a  seismic  hazard  as  is  Giles  County. 

Maximum  Giles  Countv  Earthquake 

55.  An  important  question  that  must  be  addressed  is  whether  the  largest 
possible  earthquake  has  occurred  in  Giles  County.  If  it  hasn't,  then  what  is 
the  largest  possible  earthquake  that  is  reasonable  for  this  zone?  To  answer 
this  question,  the  Waterways  Experiment  Station  engaged  Dr.  Bollinger,  a 
professor  of  seismology  at  the  Virginia  Polytechnic  Institute  and  State 
University  in  Blacksburg,  Virginia.  Dr.  Bollinger  is  a  noted  expert  on  the 
Giles  County  Earthquake  zone  and  has  published  extensively  about  the 
characteristics  of  Giles  County  seismicity.  His  report  on  the  maximum 
earthquake  potential  for  Giles  County  is  presented  in  Appendix  E. 

56.  Dr.  Bollinger  interprets  the  maximum  magnitude  earthquake  at  Giles 
County  according  to  three  different  techniques:  a)  adding  an  increment  to  the 
maximum  historical  earthquake  in  the  zone,  b)  extrapolating  the  magnitude 
recurrence  curve  for  the  zone,  and  c)  estimating  the  maximum  magnitude  as  a 
function  of  the  interpreted  fault  plane  area  that  he  has  projected  to  occur  in 
the  subsurface.  The  estimate  of  the  fault  plane  area  is  determined  primarily 
from  information  obtained  from  microearthquakes.  Bollinger's  estimation  of  M 
at  the  Giles  county  source  by  these  three  methods  are  5.9,  6.6,  and  7.0, 
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respectively  (Appendix  E) .  Bollinger's  extension  of  the  magnitude -recurrence 
curve  to  obtain  M  7.0  is  the  least  reliable  of  his  estimates  since  the 
linearity  of  the  curve  for  projection  to  large  earthquakes  is  questionable. 

An  average  of  the  remaining  two  methods  gives  M  6.3  which  is  equivalent  to 
MMI  IX.  However,  an  average  of  all  three  of  Bollinger's  methods  yields  M 
6.5  which  is  also  equivalent  to  MMI  IX.  Thus,  Bollinger's  estimates  are 
judged  to  be  consistent  with  the  MMI  IX  shown  in  Figure  12. 

Earthquake  Recurrence 

57.  A  deterministic  approach  was  use  in  this  report  to  specify 
earthquake  ground  motions.  A  deterministic  approach  is  where  a  maximum 
earthquake  is  interpreted  to  occur  at  a  fault  or  seismic  source  zone  and  the 
earthquake  is  attenuated  to  the  area  of  interest.  The  predicted  earthquake  is 
specified  for  the  seismic  region  or  zone  regardless  of  the  probability  of 
recurrence.  The  basic  assumption  is  that  the  engineered  structures  must  be 
able  to  withstand  the  predicted  intensity  of  a  maximum  credible  earthquake 
whenever  it  might  occur. 

58.  A  recurrence  relation  is  useful  for  estimating  the  general  return 
frequency  for  the  maximum  event  to  compare  to  the  operating  life  of  the 
structure.  A  recurrence  relation  is  calculated  from  the  seismic  record  for  a 
given  area  and  the  basic  Guttenburg-Richter  relationship: 

log  N  -  a  -  bM 

where  N  is  the  number  of  events  of  magnitude  M  or  greater  per  unit  of  time  and 
a  and  b  are  constants.  A  characteristic  recurrence  period  is  obtained  for  a 
given  magnitude  from  the  total  number  of  events  for  the  specified  time 
interval . 

59.  A  recurrence  relation  for  the  southeastern  United  States  and 
selected  subdivisions,  including  Giles  County,  was  developed  by  Bollinger  and 
others  (1989)  and  is  presented  in  Figures  13a  and  13b  along  with  the 
respective  recurrence  equations.  The  historical  (intensity  based)  and 
instrumental  (magnitude  based)  data  sets  were  combined  using  relations  defined 
by  Sibol,  Bollinger,  and  Birch  (1987).  The  curves  are  based  on  the  m^j^g 
magnitude  scale  (see  Appendix  A  for  description).  The  m^j^g  scale  is  considered 
generally  equivalent  to  the  m^  scale  between  nig  2  to  6.4  (Sibol,  Bollinger, 
and  Birch,  1987).  The  general  correspondence  between  itig  and  MMI  for  the 
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Eastern  United  States  Is  presented  in  Figure  14  (from  Sibol,  Bollinger,  and 
Birch,  1987). 

60.  The  mean  recurrence  for  an  MMI  VII  earthquake  for  the  southeastern 
United  States  region  is  approximately  one  every  10  years  (see  Figure  13a).  An 
MMI  VII  earthquake  is  the  threshold  where  damage  begins  to  occur  in  well  built 
engineering  structures  (see  description  of  MMI  scale.  Appendix  C) .  For  the 
Valley  and  Ridge-Blue  Ridge  provinces,  the  mean  return  period  for  an  MMI  VII 
earthquake  is  approximately  every  30  years.  For  the  Piedmont  and  the  Coastal 
Plain  provinces,  the  mean  recurrence  is  80  and  60  years,  respectively.  The 
mean  recurrence  period  for  an  MMI  VII  earthquake  in  Giles  County  is  150  years , 
for  Central  Virginia  it  is  100  years,  for  Eastern  Tennessee  it  is  80  years, 
and  at  Charleston  it  is  150  years. 

61.  The  mean  recurrence  estimated  for  a  MMI  IX  earthquake  in  Giles 
County  is  made  by  projecting  beyond  the  maximiim  historic  earthquake  recorded 
for  the  region  and  also  the  accuracy  range  of  the  data  from  which  the  curve  is 
based.  The  recurrence  estimate  for  the  maximum  earthquake  at  Giles  County 
ranges  from  310  to  4,200  years  at  a  95  percent  confidence  interval  (Bollinger 
and  others,  1989).  Dr.  Bollinger  in  his  report  in  Appendix  E  (page  E15) 
states,  "...that  the  interval  estimates,  at  a  specified  confidence  level, 
rather  than  point  estimates  are  the  preferred  manner  for  utilization  of 
magnitude  regression  results."  However,  for  purposes  of  evaluating  the 
seismic  potential  at  Gathright  Dam,  he  estimates  a  recurrence  interval  of 
approximately  1,000  years  for  the  maximum  earthquake  from  the  Giles  County 
seismic  zone. 

62.  From  the  uncertainties  in  the  recurrence  equations  and  the  general 
assumptions  that  are  made  in  the  process  of  developing  recurrence  estimates, 
the  entire  range  of  data  at  each  magnitude  interval  is  variable,  and  may 
extend  over  an  entire  log  cycle.  Because  of  this  variability  and  the 
uncertainty  with  recurrence  or  probability  estimates,  a  probabilistic  approach 
is  not  applicable  for  specifying  maximum  earthquake  ground  motions.  The 
deterministic  approach  specifies  the  maximum  credible  earthquake  regardless  of 
the  probability  of  recurrence. 
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Intensity  (MM) 


Figure  14.  Relationship  between  magnitude  and  MM  Intensity  for  eastern 
North  America;  range  in  data  is  defined  by  bar  length  and  box  plots:  mean 
equals  asterisk,  25  and  75  quartiles  equals  box  ends,  and  median  equals  center 
bar  (from  Sibol,  Bollinger,  and  Birch,  1987) 
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Felt  ’earthquakes  at  Gathright  Dam 


63.  The  southeastern  region,  with  the  exception  of  Charleston,  South 
Carolina  and  Giles  County,  Virginia,  has  low  level  earthquake  activity.  Table 
1  presents  a  list  of  MMl  VI  or  greater  earthquakes  that  were  judged  to  have 
been  felt  at  Gathright  Dam.  The  earthquake  list  is  based  on  the  catalogue  in 
Appendix  B  for  earthquakes  within  the  study  boundary  and  from  various 
published  sources  (i.e.,  Coffman,  von  Hake,  and  Stover,  1982;  Bollinger,  1975 
and  1977;  and  Bollinger  and  Hooper,  1971;  Bollinger  and  Stover,  1978;  Lessing, 
1974;  and  Stearns  and  Wilson,  1972)  for  earthquakes  centered  outside  of  the 
study  area  and  which  are  judged  to  have  been  felt  at  Gathright  Dam.  The  MMI 
at  the  source  (MMIg)  ,  the  distance  between  the  earthquake  epicenter  and 
Gathright  Dam,  and  the  attenuated  MMI  at  Gathright  Dam  (MMIg)  are  defined  for 
each  felt  earthquake  in  Table  1. 

64.  The  attenuation  procedure  selected  for  this  study  is  based  on  the 
decrease  of  intensity  with  distance  as  determined  from  curves  by  Chandra 
(1979).  The  attenuation-distance  curves  are  shown  in  Figure  15  and  the  curve 
selected  for  this  study  is  for  the  eastern  province.  The  attenuation  of  MMI 
is  determined  by  calculating  the  distance  between  the  earthquake  epicenter  and 
the  area  of  interest,  selecting  this  distance  on  the  horizontal  axis  of  the 
attenuation  curve,  and  then  deriving  the  MMI  reduction  factor.  This  reduction 
factor  is  then  subtracted  from  the  intensity  value  at  the  source  to  arrive  at 
the  interpreted  felt  intensity  at  the  site. 

65.  The  earthquakes  in  Table  1  span  approximately  200  years  and 
identify  15  events  that  were  large  enough  to  have  been  felt  at  the  damsite. 

The  vast  majority  of  felt  earthquakes  at  Gathright  Dam  have  been  less  than  MMI 
V.  There  are  five  earthquakes  in  the  list  that  are  larger  than  MMI  V,  four 
events  are  at  MMI  VI  and  one  event  is  at  MMI  VII.  There  is  only  one 
earthquake  from  these  five  that  originated  from  within  the  project  boundary. 
The  only  large  event  to  have  originated  from  inside  the  project  boundary  is 
the  1897  Giles  County  earthquake. 

66.  The  1897  Giles  County  earthquake  was  felt  at  Gathright  Dam  at  MMI 
VI  according  to  the  attenuation  distance  procedure  identified  above.  However, 
the  isoseismal  for  the  Giles  County  earthquake  shown  in  Figure  16  (from 
B'^llinger  and  Hopper,  1971)  identifies  Gathright  Dam  as  being  in  the  MMI  V 
isoseismal  area.  According  to  the  list  of  felt  intensities  for  the  towns 
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Table  1.  Felt  Earthquakes  at  Gathright  Dam 
Inside  Study  Area  Boundary  (see  Appendix  S') 

Distance  1  2 


Date 

N.  Lat 

W.  Lone 

Location 

miles 

(km) 

MM  lo 

MM  Is 

3 

May 

1897 

37.1 

80.7 

Dublin,  VA 

72 

(116) 

VII 

V 

31  May 

1897 

37.3 

80.7 

Giles  County 

60 

(97) 

VIII 

VI 

(Pearisburg,  VA) 

5 

Feb 

1898 

37.0 

80.7 

Pulaski ,  VA 

80 

(129) 

VI 

IV 

11 

Nov 

1975 

37.2 

80.9 

White  Gate,  VA 

73 

(117) 

VI 

IV 

23 

Apr 

1959 

37.4 

80.7 

Elgood,  WV 

68 

(109) 

VI 

IV 

20 

Nov 

1969 

37.4 

81.0 

Elgood,  WV 

68 

(109) 

V-VI 

III-IV 

Outside  Study  Area  Boundary 

21 

Feb 

1774 

37.3 

77.4 

Petersburg,  VA 

145 

(233) 

VII 

IV 

22 

Dec 

1875 

37.6 

77.4 

Richmond ,  VA 

140 

(225) 

VII 

IV 

23 

Dec 

1875 

37  7 

78.3 

Arvona ,  VA 

90 

(145) 

VII 

V 

10 

Oct 

1885 

37.7 

78.8 

Noorwood,  VA 

65 

(105) 

VI 

IV 

26 

Dec 

1929 

38.1 

78.5 

Charlottes¬ 

79 

(127) 

VI 

IV 

ville,  VA 

31 

Aug 

1886 

32.9 

80.0 

Charleston,  SC 

360 

(579) 

X 

II-III 

16 

Dec 

1811 

36.6 

89.6 

New  Madrid,  MO 

650 

(1046) 

XI-XII 

4 

VI 

16 

Dec 

1811 

36.6 

89.6 

New  Madrid,  MO 

650 

(1046) 

XI-XII 

VI 

23 

Jan 

1812 

36.6 

89.6 

New  Madrid,  MO 

650 

(1046) 

XI-XII 

VI 

7 

Feb 

1812 

36.6 

89.6 

New  Madrid,  MO 

650 

(1046) 

XI-XII 

4 

VII 

1.  Modified  Mercalli  Intensity  at  source  or  origin 

2.  Modified  Mercalli  Intensity  attenuated  to  site 

3.  see  Figure  18  (from  Bollinger,  1977) 

4.  see  Figure  17  (from  Stearns  and  Wilson,  1972) 
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MM  INTENSITY 
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EPICENTRAU  DISTANCE,  KM 


Figure  15.  Attenuation  of  MM  Intensity  with  distance:  A  -  Anderson,  G 
Gupta,  H-S  =  Howell  and  Schultz  (from  Chandra,  1979) 
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Figure  16.  Isoseismal  for  the  1897  Giles  County  earthquake  (from  Bollinger  and 

Hopper,  1971) 


nearest  Gathright  Dam  (from  Bollinger  and  Stover,  1978),  the  area  surrounding 
the  damsite  experienced  effects  corresponding  to  MMI  V. 

67 .  The  remaining  four  earthquakes  identified  in  Table  1  which  are 
larger  than  MMI  V  originated  from  the  New  Madrid  seismic  zone.  The  New  Madrid 
region  in  Missouri  is  the  source  for  the  four  largest  historic  earthquakes  to 
have  occurred  in  North  America.  Four  MMI  XII  earthquakes  occurred  in  1811  and 
1812  at  New  Madrid  (Street  and  Nuttli,  1984).  It  is  judged  that  the  maximum 
historic  earthquake  shaking  at  Gathright  Dam  was  experienced  during  these  New 
Madrid  events.  Stearns  and  Wilson  (1972)  identify  the  Gathright  damsite  as 
being  located  within  the  MMI  VII  isoseismal  as  shown  by  Figure  17,  a  composite 
isoseismal  for  the  New  Madrid  earthquakes  of  1811  and  1812.  Stearns  and 
Wilson  define  the  maximum  earthquake  shaking  at  the  damsite  to  have  occurred 
on  7  February  1812.  The  remaining  New  Madrid  earthquakes  are  identified  by 
Stearns  and  Wilson  and  Street  and  Nuttli  (1984)  as  being  in  the  MMI  VI  range. 

68.  The  1886  Charleston,  South  Carolina,  earthquake  produced  MMI  II  to 
III  effects  in  the  Gathright  Dam  area  as  shown  by  the  isoseismal  in  Figure  18 
(from  Bollinger,  1977).  Bollinger  and  Gilbert  (1974)  identify  the  region 
which  encompasses  the  II- III  isoseismal,  including  the  Gathright  damsite,  as 
being  in  an  earthquake  "shadow  zone,"  with  the  intensities  in  the  zone  less 
than  the  surrounding  area.  The  shadow  zone  is  part  of  an  area  that  has  been 
defined  as  having  low  heat  flow,  a  high  concentration  of  thermal  springs, 
great  crustal  thickness,  the  presence  of  Tertiary  age  volcanic  intrusives 
north  of  Gathright  Dam,  located  at  the  transition  zone  between  the  central  and 
southern  Appalachians,  and  generally  experiencing  low  level  microseismicity 
activity  (Bollinger  and  Gilbert,  1974).  The  overall  significance  of  these 
various  characteristics  is  yet  unknown.  They  suggest  that  Gathright  Dam  is 
located  in  an  area  that  experiences  low  level  seismicity  because  of  the 
underlying  geology  and  has  crustal  attenuation  characteristics  that  are 
greater  than  the  surrounding  area. 
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Figure  17.  Composite  isoseismal  for  the  1811  -  1812  New  Madrid,  Missouri, 
Earthquakes  (from  Stearns  and  Wilson,  1972) 


Figure  18.  Isoseismal  for  the  1886  Charleston,  South  Carolinea,  Earthquake 

(from  Bollinger,  1977) 
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PART  IV:  EARTHQUAKE  GROUND  MOTIONS 

Maximum  Credible  Earthquake 

69.  The  maximum  credible  earthquake  (MCE)  is  defined  as  the  largest 
possible  earthquake  that  can  be  reasonably  expected.  The  largest  earthquake 
that  is  interpreted  for  Gathright  Dam  is  an  earthquake  originating  from  the 
Giles  County  seismic  zone  (see  Figure  12).  The  maximum  earthquake  interpreted 
for  Giles  County  is  MMI  IX.  The  MCE  at  Gathright  Dam  is  the  largest  Giles 
County  earthquake  attenuated  to  the  damsite. 

70.  The  MCE  at  Giles  County  is  a  floating  earthquake  which  can  be  moved 
anywhere  within  the  source  area  shown  in  Figure  12.  However,  outside  of  this 
source  area,  the  earthquake  is  attenuated  to  the  site  of  interest  according  to 
the  distance -attenuation  relations  shown  in  Figure  15.  Gathright  Dam  is 
located  40  miles  (65  km)  from  the  edge  of  the  Giles  county  source  area.  The 
MCE  at  Giles  County,  MMI  IX,  attenuated  to  Gathright  Dam  would  produce  an 
earthquake  corresponding  to  MMI  VII,  a  reduction  of  two  intensity  levels.  An 
MMI  VII  earthquake  is  one  intensity  level  higher  than  the  seismic  region 
hosting  the  damsite  (see  Figure  12).  An  MMI  VII  earthquake  occurs  at  the 
threshold  of  damage  for  well  built  engineering  structures  (see  description  of 
MMI  in  Appendix  C) . 

71.  Ground  motions  from  earthquakes  originating  from  source  areas  other 
than  Giles  County  (see  Figure  12)  would  be  either  attenuated  with  distance  and 
would  therefore  be  less  severe  than  the  motions  from  a  maximum  Giles  County 
earthquake,  or  the  interpreted  maximum  event  for  these  source  areas  is  much 
less  than  the  estimated  maximum  for  Giles  County.  Consequently,  the  Central 
Virginia,  Northern  Virginia,  and  Southern  Appalachian  source  areas,  in 
addition  to  source  areas  that  are  more  distant,  are  not  of  concern  at 
Gathright  Dam  for  the  MCE. 

Operating  Basis  Earthquake 

72.  The  operating  basis  earthquake  (OBE)  is  an  earthquake  that  allows 
minor  damage,  but  permits  the  structure  to  remain  operational  with  small 
repairs.  It  is  an  earthquake  that  is  expected  to  occur  during  the  life  of  the 
structure.  The  life  of  the  structure  for  purposes  of  this  report  is  taken  at 
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100  years. 

73.  The  determination  of  earthquake  source  zones  for  the  OBE  at 
Gathright  Dam  is  a  difficult  consideration  in  view  of  the  low  level  of  seismic 
activity  surrounding  the  dam.  Recurrence  curves  are  presented  in  Figures  13a 
and  13b  for  the  tectonic  provinces  and  selected  seismic  hot  spots  or  source 
zones.  The  recurrence  curves  for  the  tectonic  provinces  are  conside’'ed  too 
general  for  determining  an  operating  basis  earthquake  at  Gathright  Dam.  Since 
the  dam  is  located  midway  between  several  seismic  source  zones  (see  Figure 

12) ,  it  is  more  appropriate  to  estimate  an  operating  basis  earthquake  for  each 
source  zone  and  attenuate  this  earthquake  to  the  damsite  using  the  attenuation 
and  distance  procedure  by  Chandra  (1979)  in  Figure  15. 

74.  Clearly,  the  dominant  source  area  within  62  miles  (100  km)  of 
Gathright  Dam  is  Giles  County  (see  Figure  9).  Reference  to  Figure  13b 
indicates  that  the  largest  earthquake  estimated  to  occur  at  Giles  County 
during  the  projected  100-year  operating  life  of  Gathright  Dam  is  approximately 
m^Lg  4.8,  equivalent  to  MMI  VIII  (see  Figure  14).  The  attenuated  intensity  at 
the  damsite  for  the  operating  basis  earthquake  from  a  Giles  County  source  is 
MMI  VI. 

75.  The  estimated  intensity  for  the  largest  100-year  earthquake 
occurring  in  the  Central  Virginia  seismic  zone  according  to  the  curve  in 
Figure  13b  is  m^j^^g  5.0,  equivalent  to  MMI  VIII.  This  value  is  judged  to  be 
high  considering  that  no  historic  earthquakes  larger  than  MMI  VII  have 
occurred.  Alternatively,  an  MMI  VII  earthquake  is  interpreted  for  this  zone. 
The  edge  of  the  Central  Virginia  seismic  zone  is  approximately  31  miles  (50 
km)  from  Gathright  Dam  at  its  closest  point.  As  previously  noted,  the 
determination  of  the  limits  for  hot  spots  or  seismic  source  zones  is  based  on 
the  locations  of  historic  and  microearthquakes.  The  distance  between 
Gathright  Dam  and  the  Central  Virginia  seismic  zone  represents  a  reduction  of 
the  specified  MMI  by  a  factor  of  one  according  to  the  attenuation-distance 
procedure  by  Chandra  (1979).  The  attenuated  intensity  at  the  damsite  for  the 
OBE  from  the  Central  Virginia  seismic  zone  is  MMI  VI. 

76.  The  remaining  seismic  source  area  that  has  the  potential  to  affect 
the  OBE  at  Gathright  Dam  is  the  Southern  Appalachian  seismic  zone  (see  Figure 
12).  The  Southern  Appalachian  seismic  zone  is  a  broad  belt  of  seismicity 
extending  from  Alabama  to  southwestern  Virginia.  This  zone  incorporates  Giles 
County.  The  northeastern  boundary  of  this  zone  is  located  approximately  37 


50 


miles  (60  km)  from  Gathright  Dam.  The  limits  of  the  Southern  Appalachian  zone 
correspond  approximately  to  the  Valley  and  Ridge/Blue  Ridge  Province 
identified  in  Figure  13a  by  Bollinger  and  others  (1989).  The  estimated  100- 
year  earthquake  for  the  Valley  and  Ridge/Blue  Ridge  Province  is  approximately 
nijjig  5.5,  equivalent  to  MMI  VIII.  An  MMI  VIII  earthquake  for  the  Southern 
Appalachian  zone  is  considered  high.  This  estimate  incorporates  earthquake 
data  from  both  Giles  County  and  Eastern  Tennessee,  two  areas  of  concentrated 
seismicity.  The  OBE  for  Giles  County  has  been  evaluated  separately  and  should 
be  excluded  from  this  zone  for  a  determination  of  an  OBE.  By  excluding  the 
Giles  County  zone,  the  major  center  of  seismic  activity  becomes  concentrated 
in  Eastern  Tennessee,  which  is  well  removed  from  the  damsite.  Therefore,  the 
OBE  that  is  interpreted  for  the  Southern  Appalachian  zone  is  MMI  VII,  the 
historic  maximum  for  this  broad  belt.  The  attenuated  intensity  at  the  damsite 
for  an  OBE  from  this  zone  is  MMI  V,  a  reduction  of  the  source  intensity  by  two 
intensity  levels. 

77.  It  is  concluded  that  the  maximum  OBE  at  Gathright  Dam  from  the 
various  source  areas  identified  above  is  MMI  VI.  This  value  is  equivalent  to 
the  general  maximum  defined  for  the  overall  region  surrounding  Gathright  Dam 
(see  Figure  12) . 

Field  Conditions 

78.  Ground  motions  from  an  earthquake  source  are  characterized  as  being 
either  Near  Field  or  Far  Field.  Ground  motions  are  different  for  each  field 
type.  Near  field  motions,  those  originating  near  the  earthquake  source,  are 
characterized  by  a  large  range  of  ground  motions  which  are  caused  by 
asperities  in  the  fault  plane,  complicated  reflection  and  refraction  patterns, 
and  focusing  effects  of  the  waves.  In  contrast,  the  wave  patterns  for  far 
field  motions  are  more  orderly  and  they  are  generally  muted  or  dampened. 

79.  The  limits  of  the  near  field  are  variable,  depending  on  the 
severity  of  the  earthquake.  The  relationship  between  earthquake  magnitude 
(M) ,  epicentral  intensity,  and  the  limits  of  the  near  field  are  given  in  the 
following  set  of  relations  (from  Krinitzsky  and  Chang,  1987). 
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Near  Field  Limits 


MM  Maximum 

Radius  of 

M 

Intensity.  lo 

Near  Field. 

5.0 

VI 

5 

5.5 

VII 

15 

6.0 

VIII 

25 

6.5 

IX 

35 

7.0 

X 

40 

7.5 

XI 

45 

80.  Near  field  conditions  are  specified  only  when  the  site  of  interest 
is  located  within  or  near  a  seismic  hotspot.  For  an  MMI  IX  earthquake  from 
the  Giles  County  source,  the  radius  of  the  near  field  is  22  miles  (35  km). 
Gathright  Dam  is  located  further  than  this  from  the  Giles  County  seismic 
source  area  (see  Figure  12).  A  floating  earthquake  for  the  zone  in  which 
Gathright  Dam  is  located  in  is  MM  VI .  This  zone  is  also  far  field  since  it  is 
not  a  hotspot.  Thus,  far  field  conditions  are  thus  recommended  for  the 
selection  of  ground  motions  at  Gathright  Dam. 

Recommended  Peak  Motions 


81.  The  par '.meters  for  earthquake  motions  specified  in  this  report  are 
horizontal  peak  acceleration,  horizontal  velocity,  and  duration.  Duration  is 
the  amount  of  time  in  which  the  ground  motion  is  equal  to  or  above  0.05  g 
(gravity:  1  g  =  980  cm/sec^)  .  Values  specified  are  for  free-field  motions  on 
rock  (hard  site)  at  the  surface. 

82.  The  ground  motion  parameters  of  interest  are  determined  from  the 
Krinitzsky-Chang  (1987)  intensity  curves.  The  far  field  curves  for 
acceleration,  velocity,  and  duration  are  presented  in  Figures  19,  20,  and  21. 
Values  in  the  charts  are  specified  for  the  mean,  mean  plus  one  standard 
deviation,  and  mean  plus  two  standard  deviations.  The  values  in  these  charts 
are  derived  from  a  large  world  wide  data  base  of  ground  motions  and  represent 
the  statistical  levels  for  the  spread  in  motions  for  the  different  intensity 
levels  (Krinitzsky  and  Chang,  1987). 

Maximum  Credible  Earthquake 

83.  Motions  for  the  MCE  for  Gathright  Dam  are  as  follows: 
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HOR  ACCELE 


FAR  FIELD  HARD  AND  SOFT  SITES 
ALL  MAGNITUDES 


Figure  19.  Chart  for  acceleration  (from  Krinitzsky  and  Chang,  1987) 


HOR  VELOCITY,  CM/SEC 


V  VI  VII  VIII  IX  X  XI 

MM  INTENSITY 


Figure  20.  Chart  for  velocity  (from  Krinitzsky  and  Chang,  1987) 


V  VI  VII  VIII  IX  X  XI 

MM  INTENSITY 


Figure  21.  Chart  for  duration  (from  Krinitzsky  and  Chang,  1987) 
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Giles  County  Source 

Hard  Site.  Far  Field.  MMIs  VII.  Peak  Horizontal  Motions 


Acceleration 

Velocity 

Duration 

(cm/sec^) 

(cm/sec) 

Sec.  >0.05 

Mean 

130 

9 

5 

Mean  +  S .  D . 

190 

14 

11 

84.  Where  vertical  motions  are  desired,  they  may  be  obtained  by  taking 
2/3  of  the  horizontal  values. 

Operating  Basis  Earthquake 

85.  Motions  for  an  OBE  at  Gathright  Dam  are  as  follows: 


Central  Virginia  -  Giles  Countv  Source 
Hard  Site.  Far  Field.  MMIs  VI.  Peak  Horizontal  Motions 


Acceleration 

Velocity 

Duration 

(cm/sec^) 

( cm/sec) 

Sec.  >0.05 

Mean 

80 

5 

3 

Mean  +  S.  D. 

125 

8 

6 

86.  Where  vertical  motions  are  desired,  they  may  be  obtained  by  taking 
2/3  of  the  horizontal  values. 

Recommended  Accelerograms 

87.  Three  accelerogiams  are  recommended  for  use  in  the  engineering 
analysis  of  Gathright  Dam.  The  selected  accelerograms  are  for  motions 
corresponding  to  the  mean  plus  one  standard  deviation  level.  Two  of  the 
recommended  accelerograms  are  for  the  MCE  and  the  other  accelerogram  is  for 
the  OBE.  Selected  accelerograms  are  summarized  in  Table  2  and  are  contained 
in  Appendix  F  along  with  the  associated  velocity  response  spectra,  and  the 
quadripartite  response  spectra  for  each  specified  time  history  (from  the 
California  Institute  of  Technology  Data  Base,  1975). 

88.  The  accelerograms  are  all  from  hard  sites,  a  site  in  which  the 
shear  wave  velocities  are  greater  than  1312  ft/sec  (400  meters/sec)  and  the 
underlying  geologic  horizon  is  more  than  30  ft  (9  meters)  thick.  The  scaling 
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Table 


factor  for  the  three  accelerograms  ranges  from  0.91  to  1.08.  The  scaling 
factor  is  the  ratio  between  the  recommended  peak  acceleration  and  the  peak 
acceleration  occurring  in  the  accelerogram.  Records  for  use  with  the  mean 
values  may  be  obtained  by  scaling  the  three  accelerograms  accordingly. 

Distance  from  the  source  to  the  recording  site  for  the  selected  records  ranges 
from  21  to  27  miles  (3A  to  43  km).  The  peak  motions  and  distances  are 
considered  representative  of  the  study  area. 

89.  The  records  presented  in  Table  2  are  not  the  only  records  that  may 
be  used.  Other  records  can  be  fitted  to  the  given  parameters.  The 
accelerograms  should  be  for  analogous  conditions,  such  as  size  of  earthquake, 
focal  depth  (whether  shallow  or  deep),  distance  from  source,  site  condition, 
etc.  Differences  between  peak  values  of  an  accelerogram  and  those  selected 
parameters  are  accommodated  by  changing  the  scale  of  the  accelerogram.  The 
caution  is  to  avoid  scaling  changes  that  are  greater  than  two  times  since 
larger  changes  will  affect  the  spectral  composition. 

Motions  for  Nearby  Power  Plants 


90.  The  locations  of  nuclear  and  hydroelectric  power  plants  near 
Gathright  Dam,  the  values  for  the  safe  shutdown  earthquake,  and  the  values  for 
the  operating  basis  earthquake  are  presented  in  Figure  22  (after  Nuclear  News, 
1983;  and  Dames  and  Moore,  1977).  Figure  22  identifies  three  major  power 
plants  in  Virginia.  There  are  no  such  facilities  in  West  Virginia.  The 
nearest  plant  to  Gathright  Dam  is  the  Bath  County  Pump  and  Storage  Facility. 

91.  The  safe  shutdown  earthquake  (SSE)  is  equivalent  to  the  maximum 
credible  earthquake.  Recall  that  the  OBE  is  the  maximum  earthquake  the 
structure  can  resist  and  remain  operational  without  major  damage  during  the 
design  life.  The  OBE  is  an  engineering  decision  based  on  the  cost  risk 
considerations  where  there  are  no  hazards  to  life. 

92.  Values  shown  for  peak  acceleration  for  the  SSEs  in  Figure  22  need 
not  be  directly  comparable  to  values  for  the  maximum  credible  earthquake  at 
Gathright  Dam  since  the  specification  of  values  is  dependent  on  the  types  of 
analyses  to  be  performed;  SSE  for  a  pseudostatic  analysis  would  be  a  mean 
value;  for  a  dynamic  analyses  using  an  accelerogram,  mean  plus  1  S.D.  would  be 
more  appropriate. 
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NUCLEAR  &  HYDROELECTRIC  POWER  PLANTS 

NEAR  GATHRIGHT  DAM 


ANNA 
&  4 


1  &  2 


ACCELERATION  (g)^ 


PLANT  NAME 

TYPE 

SSE  (MCE) 

OBE 

FOUNDATION 

NORTH  ANNA 
1.2,&3 

NUCLEAR 

.12 

.06 

BEDROCK 

SURRY 

1&2 

NUCLEAR 

.15 

.07 

SOIL 

BATH  CO.  PUMP 
&  STORAGE 

HYDRO 

.187 

.132 

BEDROCK 

^  SSE  -  SAFE  SHUTDOWN  EARTHQUAKE 
MCE  -  MAXIMUM  CREDIBLE  EARTHQUAKE 
QBE  -  OPERATING  BASIS  EARTHQUAKE 


Figure  22.  Locations  of  hydroelectric  and  nuclear  power  plants  and  their 
design  earthquakes  (from  Nuclear  News,  1983;  and  Dames  and  Moore,  1977) 
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Further,  the  seismic  zone  and  the  site  condition  would  introduce  other 
variations.  However,  the  motions  for  Gathright  Dam  in  Table  2  are  very  close 
to  those  obtained  independently  for  the  Bath  County  Pump  and  Storage  Facility. 
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PART  V:  CONCLUSIONS 


93.  A  seismic  zoning  was  developed  for  the  southeastern  United  States 
based  on  the  geology  and  seismic  history.  Floating  earthquakes  were  assigned 
to  each  seismic  zone  or  source  area  for  earthquakes,  since  no  active  faults 
have  been  identified  for  the  southeastern  United  States. 

94.  The  maximum  earthquake  interpreted  for  Gathright  Dam  is  from  an 
earthquake  originating  in  the  Giles  County  seismic  zone.  Gathright  Dam  is 
located  40  miles  (65  km)  from  the  Giles  County  seismic  zone.  This  zone  is 
the  location  for  the  second  largest  historic  earthquake  in  the  southeastern 
United  States.  The  maximum  credible  earthquake  at  Gathright  Dam,  attenuated 
from  the  Giles  County  seismic  zone,  is  a  far  field  earthquake  of  MMI  VII. 
Recommended  peak  horizontal  motions  for  this  earthquake  based  on  the  intensity 
curves  by  Krinitzsky  and  Chang  (1987)  are  as  follows: 


Mean 

Mean  +  S .  D . 


Maximum  Credible  Earthquake 
Hard  Site.  Far  Field.  MMIs  VII 


Acceleration 

(cm/sec^) 


Velocity  Duration 

( cm/sec)  Sec. >0.05  g 


130  9  5 

190  14  11 


95.  The  operating  basis  earthquake  interpreted  for  Gathright  Dam  is  a 
projected  100-year  earthquake  from  either  the  Giles  County  seismic  zone  or  the 
Central  Virginia  seismic  zone.  The  operating  basis  earthquake  at  Gathright 
Dam,  attenuated  from  these  source  areas,  is  a  far  field  earthquake  of  MMI  VI. 
Recommended  peak  horizontal  motions  for  this  earthquake  based  on  the  intensity 
curves  by  Krinitzsky  and  Chang  (1987)  are  as  follows: 


Mean 

Mean  +  S.  D. 


Operating  Basis  Earthquake 
Hard  Site.  Far  Field.  MMIs  VI 


Acceleration 
( cm/sec^) 


Velocity  Duration 
(cm/sec)  Sec  .^0.05  g 


80  j  3 

125  8  6 
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APPENDIX  A: 


GEOLOGY  AT  GATHRIGHT  DAM 


A1 


Detailed  information  about  the  general  geology,  foundation  geology,  dam 
design,  and  construction  characteristics  at  Gathright  Dam  are  presented  in 
various  Design  Memorandum  and/or  Construction  and  Foundation  Reports  (U.S. 

Army  Corps  of  Engineers,  1966,  1967,  1969,  1973,  1976,  1983a,  1983b,  and 
1983c) .  The  information  presented  below  and  in  other  portions  of  this  report 
is  derived  from  these  various  sources. 

Gathright  Dam  is  a  rolled  rockfill  embankment  with  an  impervious  com¬ 
pacted  earth  core  (composed  of  silty  clays  and  clayey  silts)  ,  a  double  layer 
transition  filter  (sand  and  gravel  and  quarry  spall  material)  adjacent  to  the 
core,  and  an  outer  rock  shell.  The  axis  of  the  dam  forms  an  arc  with  radius 
of  approximately  2,000  ft  (610  meters)  in  length.  The  earthen  embankment 
forming  the  dam  is  1,172  ft  (357  meters)  long  and  rises  257  ft  (78  meters) 
above  the  flood  plain  of  the  Jackson  River.  The  spillway  of  the  dam  is 
located  in  a  natural  gap  of  Morris  Hill  on  Fortney  Branch  Creek,  approximately 
2  miles  (3.2  km)  south  of  the  dam.  The  outlet  works  is  located  in  the  right 
abutment  of  the  dam  and  consists  of  a  261  ft  (80  meter)  high  reinforced 
concrete  intake  tower,  a  17.5  ft  (5.3  meter)  diameter  circular  diversion 
tunnel  extending  1181  ft  (360  meters)  through  the  right  abutment  into  a 
concrete  stilling  basin. 

The  foundations  for  the  embankment  and  abutments  are  built  on  bedrock. 
The  foundation  geology  for  the  dam  site  is  shown  in  Figure  Al  (from  U.S.  Army 
Corps  of  Engineers,  1983b).  The  rock  units  are  composed  of  Devonian  and 
Silurian  age  sediments.  The  descriptions  of  the  individual  rock  units  iden¬ 
tified  in  Figure  Al  are  presented  in  the  stratigraphic  diagram  in  Figure  A2 
(from  U.S.  Army  Corps  of  Engineers,  1976).  The  foundation  is  composed  of 
limestone,  sandstone,  and  shale.  Figure  A3  presents  a  geologic  cross  section 
along  the  centerline  of  the  dam  showing  the  distribution  of  the  different  rock 
units  (from  U.S.  Army  Corps  of  Engineers,  1983b). 

Gathright  Dam  is  located  200  to  400  ft  (approximately  60  to  90  meters) 
downstream  from  the  axis  of  the  Morris  Hill  Anticline  (see  Figure  Al) .  Con¬ 
sequently,  the  stratigraphy  at  the  damsite  dips  to  northwest  due  to  folding. 

In  general,  the  strike  and  dip  of  the  stratigraphy  in  the  reservoir  area  is 
very  irregular  because  of  the  intense  folding.  Several  different  periods  of 
deformation  are  displayed  by  the  fold  patterns. 


Foundation  geology  at  Gathright  Dam  (from  U.S.  Army  Corps  of 
Engineers,  1983b) 


SILURIAN  I  DEVONIAN 


Mlllhoro  Shale:  a  black  fissile  shale  which  weathers  light 
to  medium  gray.  Contains  concretions  locally. 


Needmore:  a  black  fissile  shale  with  concretions  and 
calcareous  siltstone  beds  locally.  Some  fossils. 


rRidqelev  or  Oriskanv;  a  reddish. buff ,  medium  grained, 

^  fossiliferous  sandstone,  contains  a  red  shale  member  in 
middle.  Sandstone  locally  conglomeratic. 

^"Becraft:  a  gray  to  blue-gray  coarse  grained  crystalline, 
fossiliferous  limestone  with  some  irregular  nodules  at 
1  base.  _  _  _ 


New  Scotland;  a  dark  gray,  argillaceous,  nodular,  chert/ 
limestone,  foss i  1 1  ferous . 


Healing  Springs:  a  buff,  fine  to  medium  grained  sandstone 
locally  calcareous  or  guar tzi tic. 

Coeywans;  blue-gray  to  pink,  coarse  grained  limestone, 
fossiliferous;  locally  an  intermediate  phase  Healing 
Springs  occurs. 

U.  Keyser:  medium  gray,  fine  to  medium  grained  limestone, 
]  wi th  undulating  bedding  and  some  chert. 

Clifton  Forge;  Light  to  medium  gray,  cross-bedded, 
calcareous  sandstone,  interbedded  with  medium  gray, 
argillaceous  limestone,  fossiliferous. 


Lower  Keyser:  Interbedded  medium  gray  argillaceous 
limestone  and  calcareous  sandstone,  fossiliferous. 


Tonoloway  Formation:  a  dark  gray,  argillaceous, 
thin  bedded,  limestone,  some  shale,  mostly  laminated; 
apparently  nonfossi 11 ferous. 


gure  A2 .  Stratigraphic  column  of  rock  units  at  Gathright  Dam  (from  U.S. 

Army  Corps  of  Engineers,  1976) 
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Figure  A3.  Geologic  cross  section  along  centerline  of  Gathright  Dam  (afte 

U.S.  Army  Corps  of  Engineers,  1983b) 


Three  major  faiilt  trends  were  identified  during  the  geological 
evaluation  of  the  damsite  and  reservoir  area  (U.S.  Army  Corps  of  Engineers, 
1976).  The  major  trends  are  based  on  the  fault  orientations.  The  faults  are 
parallel  to  the  regional  trend  (approximately  northeast-southwest), 
perpendicular  to  the  regional  trend  (approximately  east-west  to  northwest- 
southeast)  ,  and  oblique  to  the  regional  trend  (approximately  north- south) . 

The  major  faults  mapped  in  the  dam  and  reservoir  area  (see  Figure  6)  are 
generally  parallel  to  the  regional  trend.  The  mapped  faults  are  identified  as 
having  both  strike-slip  and  vertical  components  of  movement. 

There  were  no  mapped  faults  beneath  Gathright  Dam  even  though  slick-en- 
slides  were  identified  at  two  locations  beneath  the  impervious  core  foundation 
(U.S.  Army  Corps  of  Engineers,  1983b;  Plate  11-18).  At  one  location  (south 
end,  near  the  contact  between  the  Clifton  Forge  Sandstone  and  Upper  Keyser 
Limestone),  slick-en-slides  were  identified  as  being  related  to  bedding:  while 
at  the  other  location  (north  end,  near  the  contact  between  the  Clifton  Forge 
Sandstone  and  the  Lower  Keyser  Limestone) ,  there  was  no  reference  about 
stratigraphy,  only  that  slick-en-slides  are  present.  Two  northeast- southwest 
trending  faults  have  been  mapped  a  short  distance  upstream  from  the  dam.  The 
fault  nearest  to  the  dam  (see  Figures  6  and  Al)  is  identified  as  a  scissors 
fault,  where  vertical  movement  is  analogous  to  a  pair  of  scissors  with  dis¬ 
placement  in  opposite  directions  on  either  end  of  the  fault.  The  longest 
fault  identified  is  in  the  reservoir  area,  located  approximately  2000  ft  (610 
meters)  northwest  from  the  scissors  fault.  In  addition,  faults  were  mapped  on 
top  of  Hoover  Ridge  and  along  the  spillway  channel  at  the  end  of  Fortney 
Branch . 

The  jointing  of  rocks  in  the  dam  and  reservoir  area  has  influenced 
drainage  and  cave  formation  in  the  underlying  carbonate  rocks.  Four  joint 
orientations  are  present  in  the  vicinity  of  the  dam.  The  most  prominent  and 
significant  set  of  joints  are  termed  oblique  joints,  striking  between  N  65°  - 
85°  E  and  dipping  between  60°  -  85°  southeast.  Two  sets  of  dip  joints  are 
present,  striking  N  25°  -  55°  W  and  N  65°  -  85°  W  and  dipping  60°  -  90°  north¬ 
east  and  75°  -  90°  southwest,  respectively.  The  final  set  of  joints  are  iden¬ 
tified  as  bedding  and  strike  joints,  striking  N  20°  -  50°  and  dipping  50°  -  80° 
southeast  and  10°  -  70°  northwest. 

Jointing  and  solution  cavities  in  the  underlying  foundation  rocks  at  the 
dam  and  spillway  were  considered  to  be  the  major  foundation  problem  as  these 
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features  were  avenues  for  water  seepage.  Detailed  dental  work  was  performed 
in  the  foundation  of  the  dam  to  seal  these  avenues  against  possible  leakage. 
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APPENDIX  B: 


CATALOGUE  OF  HISTORIC  EARTHQUAKES 
(North  Latitude;  37.0  to  39.0,  West  Longitude:  79.0  to  81.0) 

From  Haberraann,  1989 
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APPENDIX  C:  GLOSSARY  OF  EARTHQUAKE  TERMS 


GLOSSARY 

Accelerop:ram.  The  record  Trom  an  accelerometer  presenting  acceleration  as  a 
function  of  time. 

Attenuation.  Characteristic  decrease  in  amplitude  of  the  seismic  waves  with 
distance  from  source.  Attenuation  results  from  geometric  spreading  of 
propagating  waves,  energy  absorption  and  scattering  of  waves. 

B- line .  The  slope  of  a  straight  line  indicating  frequency  of  occurrence  of 
earthquakes  versus  earthquake  magnitude. 

Bedrock .  A  general  term  for  any  hard  rock  where  it  is  not  underlain  by 
unconsolidated  materials. 

Design  Spectrum.  A  set  of  curves  used  for  design  that  shows  acceleration 
velocity,  or  displacement  (usually  absolute  acceleration,  relative  velocity, 
and  relative  displacement  of  the  vibrating  mass)  as  a  function  of  period  of 
vibration  and  damping. 

Duration  of  Strong  Ground  Motion.  The  length  of  time  during  which  ground 
motion  at  a  site  has  certain  characteristics.  Bracketed  duration  is  commonly 
the  time  interval  between  the  first  and  last  acceleration  peaks  that  are  equal 
to  or  greater  than  0.05  g.  Bracketing  may  also  be  done  at  other  levels. 
Alternatively,  duration  can  be  a  window  in  which  cycles  of  shaking  are  summed 
by  their  individual  time  intervals  between  a  specified  level  of  acceleration 
that  marks  the  beginning  and  end. 

Earthquake .  A  vibration  in  the  earth  produced  by  rupt”  e  in  the  earth's 
crust . 

1.  Maximum  Credible  Earthquake.  The  largest  earthquake  that  can  be 
reasonably  expected  to  occur. 

2.  Maximum  Probable  Earthquake.  The  worst  historic  earthquake. 
Alternatively  it  is  (a)  the  100-yoar  earthquake  or  (b)  the  -earthquake  that  by 
probabilistic  determination  of  recurrence  will  occur  during  the  life  of  the 
structure . 

3.  i^loating  Earthquake.  An  earthquake  of  a  given  size  that  can  be  moved 
anywhere  within  a  specified  area  (seismotectonic  zone). 

.  Safe  Shutdown  Earthquake.  That  earthquake  which  is  based  upon  an 
evaluation  of  the  maximum  earthquake  potential  considering  the  regional  and 
local  geology  and  seismology  and  specific  characteristics  of  local  subsurface 
material.  It  is  that  earthquake  which  produces  the  maximum  vibratory  ground 
motion  for  which  certain  structures,  systems,  and  components  are  designed  to 
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remain  functional.  These  structures,  systems,  and  components  are  those 
necessary  to  assure:  (a)  the  integrity  of  the  reactor  coolant  pressure 
boundary;  (b)  the  capability  to  shut  down  the  reactor  and  maintain  it  in  a 
safe  shutdovT?.  condition;  cr  (c)  the  capability  to  prevent  or  mitigate  the 
consequences  of  accidents  which  could  result  in  potential  offsite  exposures 
comparable  to  the  guideline  exposures  of  this  part.  (Nuclear  Regulatory 
Commission;  Title  10,  Chapter  1,  Part  100,  30  April  1975.  Same  as  Maximum 
Credible  Earthquake.) 

5.  Operating  Basis  Earthquake.  The  earthquakes  for  which  the  structure 
is  designed  to  remain  operational.  Its  selection  is  an  engineering  decision. 
Effective  Peak  Acceleration.  A  time  history  after  the  acceleration  has  been 
filtered  to  take  out  high  frequency  peaks  that  are  considered  unimportant  for 
structural  response. 

Epicenter .  The  point  on  the  earth's  surface  vertically  above  the  point  where 
the  first  earthquake  ground  motion  originates. 

Fault .  A  fracture  or  fracture  zone  in  the  earth  along  which  there  has  been 
displacement  of  the  two  sides  relative  to  one  another. 

1.  Active  Fault.  A  fault,  which  has  moved  during  the  recent  geologic 
past  (Quaternary)  and,  thus,  may  move  again.  It  may  or  may  not  generate 
earthquakes.  (U.S.  Army  Corps  of  Engineers  1983c.) 

2.  Capable  Fault.  An  active  fault  that  is  judged  capable  of  generating 
felt  earthquakes. 

Focal  Depth.  The  vertical  distance  between  the  hypocenter  or  focus  at  which 
an  earthquake  is  initiated  and  the  ground  surface. 

Focus.  The  location  in  the  earth  where  the  slip  responsible  for  an  earthquake 
was  initiated.  Also,  the  hypocenter  of  an  earthquake. 

Free  Field.  A  ground  area  in  which  earthquake  motions  are  not  influenced  by 
topography,  man-made  structures  or  other  local  effects. 

Ground  Motion.  Numerical  values  representing  vibratory  ground  motion,  such  as 
particle  acceleration,  velocity,  and  displacement,  frequency  content, 
predominant  period,  spectral  values,  intensity,  and  duration. 

Hard  Site.  A  site  in  which  shear  wave  velocities  are  greater  than  400  m/sec 
and  overlying  soft  layers  are  less  than  or  equal  to  15  m. 

Hot  Spot.  A  localized  area  where  the  seismicity  is  anomalously  high  compared 
with  a  surrounding  region. 

Intensity .  A  numerical  index  describing  the  effects  of  an  earthquake  on  man. 
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on  structures  built  by  him  and  on  the  earth's  surface.  The  number  is  rated  on 
the  basis  of  an  earthquake  intensity  scale.  The  scale  in  common  use  in  the 
U.S.  today  is  the  modified  Mercalli  (MM)  Intensity  Scale  of  1931  with  grades 
indicated  by  Roman  numerals  from  I  to  XII.  An  abridgement  of  the  scale  is  as 
follows : 

I  Not  felt  except  by  a  very  few  under  especially  favorable 
circumstances , 

II.  Felt  only  by  a  few  persons  at  rest,  especially  on  upper  floors  of 
buildings.  Delicately  suspended  objects  may  swing. 

III.  Felt  quite  noticeable  indoors,  especially  on  upper  floors  of 
buildings,  but  many  people  may  not  recognize  it  as  an  earthquake.  Standing 
motor  cars  may  rock  slightly.  Vibration  like  passing  of  truck.  Duration  can 
be  estimated. 

IV.  During  the  day  felt  indoors  by  many,  outdoors  by  few.  At  night 
some  awakened.  Dishes,  windows,  doors  disturbed;  walls  make  cracking  sound. 
Sensation  like  heavy  truck  striking  building.  Standing  motor  cars  rocked 
noticeably . 

V.  Felt  by  nearly  everyone;  many  awakened.  Some  dishes,  windows, 
etc.,  broken;  a  few  instances  of  cracked  plaster;  unstable  objects  overturned. 
Disturbance  of  trees,  poles  and  other  tall  objects  sometimes  noticed. 

Pendulum  clocks  may  stop. 

VI.  Felt  by  all;  many  frightened  and  run  outdoors.  Some  heavy 
furniture  moved;  a  few  instances  of  fallen  plaster  or  damaged  chimneys. 

Damage  slight. 

VII.  Everybody  runs  outdoors.  Damage  negligible  in  buildings  of  good 
design  and  construction;  slight  to  moderate  in  well-built  ordinary  structures; 
considerable  in  poorly  built  or  badly  designed  structures;  some  chimneys 
broken.  Noticed  by  persons  driving  motor  cars. 

VIII.  Damage  slight  in  specially  designed  structures;  considerable  in 
ordinary  substantial  buildings  with  partial  collapse;  great  in  poorly  built 
structures.  Panel  walls  thrown  out  of  frame  structures.  Fall  of  chimneys, 
factory  stacks,  columns,  monuments,  walls.  Heavy  furniture  overturned.  Sand 
and  mud  ejected  in  small  amounts.  Changes  in  well  water.  Persons  driving 
motor  cars  disturbed. 

IX.  Damage  considerable  in  specially  designed  structures;  well 
designed  frame  structures  thrown  out  of  plumb;  damage  great  in  substantial 
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buildings,  with  partial  collapse.  Buildings  shifted  off  foundations.  Ground 
cracked  conspicuously.  Underground  pipes  broken. 

X.  Some  well-built  wooden  structures  destroyed;  most  masonry  and 
frame  structures  destroyed  with  foundations;  ground  badly  cracked.  Rails 
bent.  Landslides  considerable  from  river  b.-inks  end  steep  slopes.  Shifted 
sand  and  mud.  Water  splashed  over  banks. 

XI.  Few  structures  remain  standing.  Unreinforced  masonry  structures 
are  nearly  totally  destroyed.  Bridges  destroyed.  Broad  fissures  in  ground. 
Underground  pipe  lines  completely  out  of  service.  Earth  slumps  and  land  slips 
in  soft  ground.  Rails  bent  greatly. 

XII.  Damage  total.  Waves  seen  on  ground  surfaces.  Lines  of  sight  and 
level  distorted.  Objects  thrown  upward  into  the  air. 

Liquefaction.  The  sudden,  total  loss  of  shear  strength  in  a  soil  as  the 
result  of  excess  pore  water  pressure.  The  result  is  a  temporary 
transformation  of  unconsolidated  materials  into  a  fluid. 

Magnitude .  A  measure  of  the  size  of  an  earthquake  related  to  the  strain 
energy.  It  is  based  upon  the  displacement  amplitude  and  period  of  the  seismic 
waves  and  the  distance  from  the  earthquake  epicenter. 

1.  Body  Wave  Magnitude  Cm^')  .  The  mj,  magnitude  is  measured  as  the  common 
logarithm  of  the  maximum  displacement  amplitude  (microns)  of  the  P-wave  with 
period  near  one  second.  Developed  to  measure  the  magnitude  of  deep  focus 
earthquakes,  which  do  not  ordinarily  set  up  detectable  surface  waves  with  long 
periods.  Magnitudes  can  be  assigned  from  any  suitable  instrument  whose  con¬ 
stants  are  known.  The  body  waves  can  be  measured  from  either  the  first  few 
cycles  of  the  compression  waves  (mj,)  or  the  1  second  period  shear  waves 
(niblg)  • 

2.  Local  Magnitude  (M,^)  .  The  magnitude  of  an  earthquake  measured  as  the 
common  logarithm  of  the  displacement  amplitude,  in  microns,  of  a  standard 
Wood-Anderson  seismograph  located  on  firm  ground  62  miles  (100  km)  from  the 
epicenter  and  having  a  magnification  of  2,800,  a  natural  period  0.8  second, 
and  a  damping  coefficient  of  80  percent.  Empirical  charts  and  tables  are 
available  to  correct  to  an  epicentral  distance  of  62  miles  (100  km),  for  other 
types  of  seismographs  and  for  various  conditions  of  the  ground.  The 
correction  charts  are  suitable  up  to  epicentral  distances  of  373  miles  (600 
km)  in  southern  California  and  the  definition  itself  applies  strictly  only  to 
earthquakes  having  focal  depths  smaller  than  about  19  miles  (30  km).  The 
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correction  charts  are  suitable  up  to  epicentral  distances  of  about  373  miles 
(600  km) .  These  correction  charts  are  site  dependent  and  have  to  be  developed 
for  each  recording  site. 

3.  Surface  Wave  Magnitude  (M^).  This  magnitude  is  measured  as  the 
common  logarithm  of  the  resultant  of  the  maximum  mutually  perpendicular 
horizontal  displacement  amplitudes,  in  microns,  of  the  20-second  period 
surface  waves.  The  scale  was  developed  to  measure  the  magnitude  of  shallow 
focus  earthquakes  at  relatively  long  distances.  Magnitudes  can  be  assigned 
from  any  suitable  instrument  whose  constants  are  known. 

4.  Richter  Magnitude  (M) .  Richter  magnitude  is  nonspecified  but  is  usu¬ 
ally  M^^  up  to  6.5  and  Mg  for  greater  than  6.5. 

5.  Seismic  Movement  (M^)  .  Seismic  moment  is  an  indirect  measure  of 
earthquake  energy. 

M„  =  G  A  D 

where 

G  =  rigidity  modulus 

A  =  area  of  fault  movement 

D  =  average  static  displacement 
The  values  are  in  dyne  centimeters. 

6.  Seismic  Moment  Scale  (M.. )  .  Expresses  magnitude  based  on  the  concept 
of  seismic  moment: 


M^  =  2/3  log  M„  -  10.7 


7.  Comparison  of  Magnitude  Scales.  Table  7-1  presents  a  comparison  o<^ 
values  for  m^,  M^,  M,  log  M^,  M,„  and  Mg. 


Table  7-1. 

Comparison 

between 

mt,.  Ml,  M,  log  M„, 

M,,,  and 

Mg  scales . 

n't 

Ml 

M 

LogM^j  ( dyne  -  cm ) 

M 

Ms 

Bodv-Wave 

Local 

Richter 

Seismic  Moment 

Moment 

Surface -Wave 

5.0 

5.4 

5.4 

24.2 

5.4 

5.0 

5.5 

5.9 

5.9 

25.0 

6.0 

5.8 

6.0 

6.4 

6.7 

26.1 

6.7 

6.7 

6.5 

6.9 

7 . 5 

27.3 

7.5 

7  .  5 

7.0 

7  .  5 

8.3 

28.6 

8.4 

8.3 

Particle  Acceleration.  The  time  rate  of  change  of  particle  velocity. 
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Particle  Displacement.  The  difference  between  the  initial  position  of  a 
particle  and  any  later  temporary  position  during  shaking. 

Particle  Velocity.  The  time  rate  of  change  of  particle  displacement. 

Response  Spectrum.  The  maximum  values  of  acceleration,  velocity,  and/or 
displacement  of  an  infinite  series  of  single-degree-of- freedom  systems,  each 
characterized  by  its  natural  period,  subjected  to  a  time  history  of  earthquake 
ground  motion.  The  spectrum  of  maximum  response  values  is  expressed  as  a 
function  of  natural  period  for  a  given  damping.  The  response  spectrum 
acceleration,  velocity,  and  displacement  values  may  be  calculated  from  each 
other  by  assuming  that  the  motions  are  harmonic.  When  calculated  in  this 
manner  these  are  sometimes  referred  to  as  pseudo-acceleration,  pseudo¬ 
velocity,  or  pseudo-displacement  response  spectrum  values. 

Saturation.  Where  those  measures  of  earthquake  motions  (acceleration, 
velocity,  magnitude,  etc.)  do  not  increase  though  the  earthquakes  generating 
them  may  become  larger. 

Scaling .  An  adjustment  to  an  earthquake  time  history  or  response  spectrum 
where  the  amplitude  of  acceleration,  velocity,  and/or  displacement  is 
increased  or  decreased,  usually  without  change  to  the  frequency  content  of  the 
ground  motion. 

Seismic  Hazard.  The  physical  effects  of  an  earthquake. 

Seismic  Risk.  The  probability  that  an  earthquake  of  or  exceeding  a  given  size 
will  occur  during  a  given  time  interval  in  a  selected  area. 

Seismic  Zone.  A  geographic  area  characterized  by  a  combination  of  geology  and 
seismic  history  in  which  a  given  earthquake  may  occur  anywhere. 

Soft  Site.  A  site  in  which  shear  wave  velocities  are  less  than  400  m/sec  in  a 
surface  layer  16  or  moi.e  m  thick. 
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APPENDIX  D: 


IN3TRUMENTALLY  LOCATED  EARTHQUAKES  IN  VIRGINIA 
(From  Bollinger  and  others,  1986) 
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EXECUTIVE  SUMMARY 


The  maximum  magnitude  earthquake  expected  from  the  Giles  county, 
Virginia  seismic  zone  is  estimated.  Of  the  various  techniques  employed  to 
obtain  estimates  for  such  maximum  seismic  events,  three  are  applicable 
to  the  data  bases  available  for  eastern  United  States  seismicity  ;  (1) 
Adding  an  increment  to  the  maximum  historical  earthquake  in  the  zone,  (2) 
Extrapolation  of  the  magnitude  recurrence  curve  for  the  zone,  and,  (3) 
Magnitudes  based  on  estimates  of  the  fault  area  of  the  zone.  Each  of  these 
techniques  has  associated  uncertainties  both  in  their  applicability  to  the 
zone  under  consideration  as  well  as  in  the  determination  of  the  key 
parameters  involved.  The  process  of  maximum  magnitude  estimation  is 
intrinsically  subjective  and  depends  directly  on  the  experience  and 
judgment  of  the  analyst. 

Application  of  the  above  three  techniques  to  the  Giles  county, 
Virginia,  seismic  zone  leads  to  the  following  results  ; 

Ms  =  6.9  from  adding  a  1.0  increment  to  the  maximum  historical 
earthquake  known  to  have  occurred  in  the  zone  (May  31,  1897  ;  MMI  =  VIII  : 
mb  =  5.8,  Ms  =  5.9), 

Ms  =  6.95  from  extension  of  the  magnitude-recurrence  curve,  and 

Ms  =  6.57  from  the  average  of  six  estimates  for  the  fault  zone  area 
ranging  from  112  sq  km  (Ms  =  6.34)  to  300  sq  km  (Ms  =  6.76). 

For  a  single  estimate  of  maximum  magnitude  the  average  of  the  above 
three  values,  rounded  to  the  nearest  one-tenth,  should  be  used.  That  value 
is  : 

Ms  =  6.8  or  equivalently,  mb  =  6.3  . 

For  multiple  estimates,  the  two  extreme  values  can  be  utilized. 


The  Definition  of  Maximum  Magnitude  Earthquake 

Possible  wavs  to  define 

Some  of  the  ways  in  which  maximum  magnitude  have  been  defined  are: 
(1)  The  largest  possible  earthquake  that  can  occur  given  the  current 
physical  conditions  (no  change  in  the  future)  of  the  source  area,  or,  (2)  The 
largest  possible  earthquake  to  occur  with  a  specified  probability  during  a 
specified  exposure  time,  or,  (3)  The  largest  earthquake  likely  to  occur  in  a 
reasonable  amount  of  time  (life  of  facility  involved  ?).  Note  that  (3)  is  a 
qualitative  form  of  (2). 

Possible  synonyms 

Maximum  Possible  Earthquake  ;  Maximum  Credible  Earthquake.  The 
point  here  that  these  different  terms  can  have  various  meanings  to 
different  individuals.  The  expression  'maximum  magnitude  earthquake'  will 
be  used  herein  and  will  be  defined  subsequently. 

Problems 

Given  the  long  recurrence  intervals  for  the  larger  intraplate 
earthquakes  and  the  short  historical  record,  there  is^  the  possibility  that 
the  "maximum"  earthquake  has  not  been  recorded  in  a  given  zone.  Poisson 
statistics  would  indicate  that  there  is  a  63%  probability  of  an  earthquake 
catalog  containing  an  earthquake  with  recurrence  interval  equal  to  or 
greater  than  the  length  of  the  catalog.  This  suggests  that,  in  more  than 
one-third  of  the  earthquake  catalogs,  there  is  likely  to  be  an  apparent 
deficiency  of  large  shocks  (Chinnery,1979).  Another  way  to  utilize  Poisson 
statistics  in  this  instance  :  In  the  southeastern  U.S.,  the  earthquake  record 
is  about  250  years  long.  There  is  a  1  in  5  chance  (22%)  that  it  contains  a 
shock  with  a  recurrence  interval  of  1000  years  which  has  been  suggested 
as  a  candidate  definition  for  maximum  magnitude. 
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Chinnery  (1979)  has  demonstrated  that  there  is  no  proof  for  an 
"absolute"  upper  bound  to  seismic  moment,  and  hence  earthquake  size,  on  a 
global  scale  even  though  there  are  physical  arguments  that  such  an  upper 
bound  must  exist.  Thus,  he  points  out  that  there  is  no  unequivocal  way  to 
know  with  certainty  if  you  have  the  maximum  earthquake  in  a  given 
catalog  or  not. 

Definition  for  Maximum  Magnitude  for  this  Study. 

The  purpose  of  this  study  is  to  estimate  the  maximum  magnitude 
earthquake  for  the  Giles  county,  Virginia,  seismic  zone.  Accordingly,  we 
adopt  a  definition  similar  to  (2)  and  (3)  above.  The  "specified  probability" 
will  be  set  at  0.001,  without  consideration  of  exposure  time.  The  question 
of  whether  or  not  such  a  seismic  event  can  occur  elsewhere  in  the  area  is 
not  addressed. 


The  Estimation  of  Maximum  Magnitude 


Use  of  the  Historical  Record  of  Earthquakes 

Recurrence  relationships  assume  that  the  past  record  of  small  and 
large  earthquakes  is  representative  of  future  seismic  activity  for  as  long 
as  is  necessary.  However,  very  long  catalogs  from  seismically  active 
interplate  or  plate  marginal  areas,  e.g.,  the  Middle  East,  China,  and  Japan, 
show  long  term  changes  in  seismicity  on  time  scales  of  100's  of  years. 
Whether  or  not  such  secular  variations  are  also  appropriate  for  intraplate 
settings  is  not  known  for  certain.  Thus,  the  possibility  exists  that  the 
largest  earthquakes  for  a  given  seismic  zone  may  be  associated  with  a 
level  of  seismicity  that  is  very  different  from  the  recent  record  there  of 
smaller  shocks. 

The  recurrence  relation  (Log  N  versus  M)  must  be  related  to  a 
maximum  magnitude  in  some  manner.  If  a  physical  limit  or  a 
"characteristic"  earthquake  does  exist  for  a  given  zone,  then  the 
recurrence  curve  will  have  to  cut-off  abruptly  or  bend  rapidly  in  some 
manner  so  as  to  be  parallel  to  the  ordinal  axis  at  that  magnitude.  This 


factor  impacts  the  simple  extension  of  a  magnitude  recurrence  curve  to 
larger  magnitudes  in  the  maximum  magnitude  estimation  process. 


In  terms  of  the  maximum  historical  earthquake  for  the  eastern  U.S. 
host  region,  there  were,  e.g.,  the  Ms  >  8  shocks  at  New  Madrid,  Missouri, 
and  the  Ms  >  7  shock  at  Charleston,  South  Carolina.  Thus,  at  least  two 
locations  in  the  region  have  exhibited  moderate  to  large  earthquakes.  New 
Madrid  and  Charleston  are  also  the  only  seismic  areas  east  of  the  Rocky 
Mountains  that  have  paleoseismic  evidence  for  pre-historical  occurrences 
of  larger  shocks.  There  were  3  major  earthquakes  in  the  past  2000  yearc 
or  less  at  New  Madrid  (Russ,  1979)  and  3  moderate  or  larger  shocks  in  the 
past  7200  years  at  Charleston  (Obermier  and  others,  1987).  When  the 
recurrence  curves  for  those  areas  (excluding  the  largest  historical  shocks) 
were  projected  to  repeat  times  of  some  600  years  in  Missouri  and  1000 
years  in  South  Carolina,  the  magnitudes  indicated  were  in  good  agreement 
with  the  estimated  magnitudes  for  the  largest  historical  earthquakes 
(Nuttli,  1981). 

The  agreement  of  historical  and  pre-historical  data  in  Missouri  and 
South  Carolina  is  very  important  to  the  estimation  of  maximum  magnitude 
earthquakes  in  the  region  as  both  shocks  are  large  enough  to  be  reasonable 
candidates  for  the  maxima  in  their  respective  zones.  Under  the  assumption 
that  such  is  the  case  we  have  :  (1)  Different  seismic  zones  in  the  eastern 
U.S.  can  have  different  maximum  magnitude  earthquakes,  i.e.,  some  zones 
have  smaller  maximum  magnitudes  than  other  zones,  and,  (2)  The  rate  of 
strain  accumulation,  amount  of  fault  surface,  and  the  friction  on  the  fault 
surfaces  are  different  for  different  source  volumes. 

Magnitude  Recurrence  Relations  for  the  Giles  Countv  Seismic  Zone  fGCSZ) 

Bollinger  and  others  (1989a)  have  recently  completed  an  extensive 
study  of  frequency  of  earthquake  occurrence  in  the  southeastern  U.  S.  That 
study  included  investigation  of  the  Giles  county,  Virginia,  seismic  zone. 
Their  results  will  be  utilized  in  this  study. 

Specification  of  the  Area  of  the  GCSZ. 

In  some  seismic  hazard  studies,  it  is  necessary  to  normalize  for  the 


area  (volume)  being  considered.  Otherwise,  there  would  be  no  limit  (other 
than  global)  to  how  large  a  magnitude  could  be  estimated  as  larger  and 
larger  source  regions  are  considered  (Nuttli,1981).  However,  the  GCSZ  is 
small  enough  (7,854  sq  km  (Davison,  1988  ;  Bollinger  and  others,  1989a,b) 
that  no  normalization  is  required  for  the  task  at  hand. 

The  definition  of  the  actual  boundaries  of  the  seismic  source  zone  is 
not  without  its  own  uncertainties,  because  most  active  sources  tend  to 
display  a  'halo'  of  surrounding  seismic  activity.  That  activity  is  generally 
assumed  to  be  due  tn  peripheral  stress  perturbations  induced  by  the  zone 
proper.  Such  halos  blur  the  actual  boundaries  of  the  principal  zone, 
especially  given  errors  in  hypocentral  locations.  There  is  also  the  very 
real  question  as  to  whether  or  not  the  halo  activity  should  be  considered 
an  integral  part  of  the  zone. 

The  GCSZ  has  the  distinction  of  being  the  site  for  the  second  largest 
earthquake  known  to  have  occurred  in  the  southeastern  U.S.  Its 
meizoseismal  intensity  was  MMI  VIII  and  Nuttli  and  others  (1979,  1989) 
have  estimated  magnitudes  of  mb  =  5.8  and  Ms  =  5.9.  Its  small 

meizoseismal  area  indicates  an  epicenter  near  the  county  seat  of 
Pearisburg  (Bollinger  and  Hopper,  1971).  The  spatial  distribution  of  the 
historical  seismicity  (Bollinger,  1973a,b)  shows  the  zone  to  be  relatively 
isolated,  but  not  sharply  defined.  The  results  from  a  decade  of  monitoring 
by  a  seismic  network  sited  to  study  the  zone  have  corroborated  the 

principal  historical  results  (Bollinger  and  others,  1986)  that  there  is  an 
area  of  isolated  seismicity  in  the  Giles  county  locale,  but  its  spatial 
configuration  is  not  simple  (Bollinger  and  others,  1989b). 

Fortunately,  eight  well  constrained  sets  of  focal  mechanism 

solutions,  based  on  both  P-wave  polarity  and  S/P  wave  amplitude  ratios, 
have  been  developed  for  Giles  county  earthquakes  (Munsey  and 

Bollinger, 1 985  ;  Davison, 1 988).  Those  focal  mechanisms  were  used  by 
Davison  (1988)  to  estimate  the  regional  in-situ  stress  as  being 
compressive  and  northeasterly  trending.  Given  that  estimate,  Davison 
(1988)  was  then  able  to  select  a  preferred  fault  plane  from  each  pair  of 
nodal  planes  on  the  basis  of  compatibility  between  the  direction  of  the 
slip  on  each  nodal  plane,  as  indicated  by  the  focal  mechanism,  and  the 
direction  of  slip  expected  from  the  regional  stresses.  Given  the 


orthogonal  relationship  between  the  nodal  plane  pairs  for  each  earthquake, 
slip  compatibility  with  the  regional  stresses  is  an  effective  criterion. 
Davison  (1988)  determined  the  average  of  the  preferred  nodal 
plane  strikes  to  be  N25°E.  It  is  important  to  note  that  this  average 
strike  estimate  is  based  solely  on  focal  mechanism  data  and  is 
independent  of  any  direct  interpretation  of  epicentral  patterns,  the 
technique  usually  employed  to  identify  earthquake  fault  zones. 

The  entire  earthquake  catalog  for  the  GCSZ  is  shown  in  Figure  1  and 
listed  in  Appendix  A.  When  the  very  poorly  constrained  epicenters  for 
historical  and  recent  shocks  are  deleted  so  as  to  leave  only  those  whose 
hypocenters  are  known  within  ±  10  km,  the  pattern  that  remains  is  shown 
by  Figure  2.  Utilizing  this  data  set  (listed  in  Appendix  B)  in  conjunction 
with  the  average  strike  of  N25°E  will  allow  the  area  (volume)  of  the  GCSZ 
to  be  estimated.  Figures  3  through  5  show  the  GCSZ  definition  based  on  a 
±10  km  width  on  either  side  of  the  N25°E  trending  line  through  the 
hypocentral  lineation.  That  definition  provides  a  geological  interpretation 
of  the  seismic  observational  results. 


Use  of  Fault  Plane  Area  -  Magnitude  Relationships 

The  data  bases  for  such  relationships  are  almost  entirely  from 
interplate  and  plate-marginal  regions  that  are  very  active  seismicaily 
(high  strain  rates)  and  often  exhibit  surface  faulting  associated  with  the 
causal  faults.  The  applicability  of  such  results  to  a  low  activity, 
intraplate  region  containing  only  buried  causal  faults,  some  at  relatively 
large  depths,  is  questionable.  The  data  bases  themselves  are  not  without 
some  questions  as  to  their  adequacy  and  quality.  However,  the  basic 
physics  of  the  seismogenic  process  is  contained  in  the  spatial  fault 
parameters  and  some  of  the  fault  areas  were  estimated  with  the  help  of 
aftershock  surveys.  They  can,  therefore,  be  used  effectively  as  part  of  the 
maximum  magnitude  estimation  procedure. 

The  physical  theory  of  the  earthquake  process  indicates  that 
earthquake  magnitude  should  be  more  strongly  correlated  with  the 
logarithm  of  the  fault  area  than  with  the  logarithm  of  the  fault  length 
alone.  Wyss  (1979,  1980  ;  Bonilla,  1980)  and  Singh  and  others  (1980)  have 
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Figure  1.  Seismicity  Maps  for  the  Giles  county,  Virginia,  seismic  zone  -  1876  through  1988. 
Epicenters  indicated  by  octagon  symbols.  Circular  definition  of  zone  (radius  =  50  km', 
according  to  Davison,  1988.  N  =  number  of  epicenters  plotted.  Upper  figure  :  Epicenters  only 
Lower  figure  :  Same  epicenters  with  horizontal  error  bars.  The  large  error  bars  are  for 
historical  shocks  for  which  instrumental  control  was  lacking  or  sparse.  The  small  error  bars 
are  the  result  of  monitoring  by  a  local  network  of  seismographs  (Bollinger  et  al.  1986). 
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Figure  2.  Seismicity  Map  for  the  Giles  county,  Virginia,  seismic  zone  -  1876  through  1988 
showing  only  those  earthquakes  for  which  the  epicenters  and  focal  depths  are  known  within 
±  10  km.  Symbols  and  format  are  the  same  as  in  Figure  1. 
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Figure  3.  Tabular  definition  ‘or  the  Giles  county.  Virginia,  seismic  zone.  Epicenter  and  error 
bar  ^mbols  same  as  in  Figure  1.  Upper  figure  :  Dashed  box  at  10  km  about  a  line  trending 
N25‘^E  encloses  the  zone  proper  ;  Lower  figure  ;  Epicenters  within  the  dashed  box  shown  in 
the  upper  figure  are  projected  into  the  vertical  plane  trending  N250e  shown  here.  Both 
horizontal  and  vertical  error  bars  shown  for  each  focus.  Profile  distance  measured  from 
southwest  to  northeast. 
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Figure  4,  Tabular  defini'ion  for  the  Giles  county,  Virc'nia,  seismic  zone.  Epicenter  and  emor 
bar  ^mpols  same  as  in  Figure  1.  Upper  figure  :  Cashed  box  at  ±10  km  aboui  a  line  trending 
N25<JE  encloses  the  zone  proper  :  Lower  figure  :  Epicenters  within  the  dashed  box  shown  in 
the  upper  figure  projected  into  the  plane  perpendicular  to  the  N25^E  trend  shgwn  here.  Both 
horizontal  and  vertical  error  bars  shown  for  each  focus.  Profile  distance  measured  from 
northwest  to  southeast. 
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Figure  5.  Seismicity  map  for  the  Giles  county,  Virginia,  seismic  zone  -  1876  through  1988, 
showing  all  epicenters  (octagon  symbols)  and  the  circular  and  rectangular  definitions  of  the 
zone.  The  circular  definition  is  a  general  one  and  includes  the  off-zone  'halo'  events.  The 
rectangular  definition  is  for  the  zone  proper  and  includes  only  those  earthquakes  thought  to 
have  originated  within  the  principal  seismogenic  structure. 


developed  regressions  of  magnitude  (M  =  mostly  Ms  with  some  Ml  and  Mw  ; 
M  >  5.6)  with  fault  plane  area  (A,  sq  km)  .  Bonilla  and  others  (1984)  prefer 
the  use  of  fault  length  (L)  and/or  displacement,  but  they  also  presented 
magnitude-fault  area  regressions.  These  equations  are  as  follows  : 

Wyss  (1979)  M  =  4.15  +1.00  Log  A 

Singh  et  al.  (1980)  M  =  4.53  +  0.89  Log  A 

Bonilla  et  al.  (1984)  M  =  4.36  +  1.035  Log  A 

These  expressions  are  similar  to  each  other  (they  used  similar  data  bases)^ 
and  there  is  no  obvious  basis  to  prefer  any  one  of  them  over  the  others  and, 
thus,  the  average  of  their  results  will  be  employed  herein.  Furthermore, 
the  M  will  be  interpreted  as  Ms,  because  of  the  preponderance  of  that 
earthquake  size  measure  in  the  input  data  bases.  Bonilla  and  others  (1984) 
magnitude  -  fault  length  expression,  M  =  6.02  +  0.729  Log  L,  will  be 
utilized  only  in  a  comparative  manner  as  it  relates  primarily  to  surface 
ruptures. 

Estimation  Procedures  and  Results  for  the  Giles  County, 

Virginia  Seismic  Zone 

The  problems  discussed  in  the  preceding  section  must  be  dealt  with 
by  the  analyst  on  a  case  by  case  basis.  The  decision  must  be  made  as  to 
which  techniques  are  applicable  to  the  area  being  studied.  Some  problems, 
e.g.,  the  fact  that  a  global  maximum  earthquake  has  not  yet  been 
documented,  can  only  be  recognized.  For  other  problems,  e.g.,  the  use  of 
interplate  magnitude-fault  area  relations  in  intraplate  environments,  it  is 
necessary  for  the  analyst  to  present  the  judgments  and  reasoning  utilized 
to  justify  their  use  or  non-use  on  the  study  area  being  considered. 

Applications  of  Historical  Seismicity 

Increment  to  the  Historical  Maximum  Earthquake 

In  practice,  0.5  or  1.0  magnitude  units  have  been  sometimes  been 
added  to  the  largest  historical  earthquake  as  an  estimate  of  the  maximum 
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shock  for  a  zone.  This  is  a  subjective  procedure  that  depends  completely 
on  the  judgment  of  the  analyst.  The  only  quantitative  aspect  of  this 
procedure  is  the  fact  that  such  an  addition  actually  implies  an  assumed 
lengthening  of  the  historical  record.  Thus,  for  a  b-value  of  -1,  a  0.5 
addition  implies  a  3.2  times  lengthening  of  the  historical  record,  while  a 
1.0  increment  implies  a  factor  of  10  times.  Recognition  of  the  actual 
amount  of  time  extension  should  be  made  in  each  particular  case. 

Results  for  the  GCSZ 

For  the  GCSZ  the  b-value  is  0.64  (Bollinger  and  others,  1989).  Thus,  an 
increment  of  0.5  magnitude  units  implies  a  factor  of  2.09X  and  a  1.0 
increment  a  4.37X  factor.  The  earthquake  catalog  for  the  GCSZ  is  215  years 
and  those  factors  imply  extension  time  intervals  of  450  years  and  940 
years,  respectively.  As  noted  previously,  the  maximum  shock  for  the  zone 
was  a  Ms  =  5.9.  An  increment  of  1.0  Ms  units  is  selected  as  both  a 
conservative  measure  and  as  one  that  is  compatible  with  the  definition  of 
maximum  magnitude  adopted  herein.  Thus,  the  estimation  for  this 
procedure  is  Ms  =  6.9  which  implies  an  extension  time  of  940  years. 

Extrapolation  of  the  Recurrence  Curve 

This  procedure  is  similar  to  the  preceding  one  except  that  the 
objective  is  a  given  recurrence  interval,  rather  than  a  given  magnitude 
increment.  That  is,  the  recurrence  curve  extrapolation  results  are  directly 
dependent  on  the  specific  intercept  (  a  )  and  slope  (  -b  )  values  of  the 
curve  being  extended.  The  magnitude  increment  addition  is  completely 
independent  of  the  a  and  b  values. 

These  extrapolations  are  usually  linear,  but  it  is  well  documented  in 
the  western  U.S.  that  the  difference  between  'background  seismicity'  and 
large,  'characteristic  earthquakes’  is  nonlinear.  However,  linear 
extrapolation  is  the  most  conservative  with  respect  to  the  various 
truncated  or  exponential  fall-off  terminations  proposed  for  log  N  versus  M 
curves  and  will  be  employed  herein. 

Nuttli  (1981)  recommends  use  of  the  magnitude  associated  with  a 
1000  year  recurrence  interval  (annual  probability  =  0.001)  for  seismic 
source  zones  (normalized  to  30,000  (or  less)  sq  km  or  100,000  sq  km)  as 
an  estimate  of  the  maximum  magnitude  for  eastern  U.S.  source  zones.  His 
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recommendation  is  based  on  analyses  of  the  Mew  Madrid,  Missouri  and 
Charleston,  South  Carolina  zones.  In  both  of  those  zones,  he  deleted  the 
largest  historical  earthquakes,  determined  a  magnitude-recurrence 
relation  for  the  remaining  catalog,  and  then  extrapolated  the  resulting  Log 
N  versus  M  curve  to  a  recurrence  interval  of  1000  years.  The  magnitudes 
associated  with  the  1000  year  intervals  were  in  good  agreement  with 
those  for  the  largest  historical  events  that  had  been  deleted.  A  problem 
with  Nuttli's  (1981)  approach  is  that  the  recurrence  relationships  must  be 
normalized  to  some  arbitrary  area  to  yield  consistent  results. 
Furthermore,  its  applicability  to  seismic  zones  with  very  small  areas, 
such  as  the  GCSZ,  has  not  been  demonstrated.  Acknowledging  those 
problems,  we  choose  the  1000  year  earthquake  for  the  zone  as  a 
reasonable  estimate  of  the  GCSZ  maximum  magnitude  shock. 

Nuttli  (1981)  also  noted  that,  "East  of  the  Appalachians,  the 
earthquake  source  zones  are  not  so  readily  delineated,  so  it  is  difficult  to 
assign  maximum  magnitude  earthquakes  to  that  part  of  the  country."  We 
agree  with  that  assessment. 

Results  for  the  GCSZ 

Bollinger  and  others  (1989a)  have  determined  the  recurrence 
relationship  for  the  GCSZ  as. 

Log  Nc  =  1.065  -  0.64  mb(Lg). 

That  equation  yields  a  mb(Lg)  =  6.35  (ms  =  6.95)  for  a  recurrence  interval 
of  1000  years.  Bollinger  and  others  (1989)  note  that  interval  estimates,  at 
a  specified  confidence  level,  rather  than  point  estimates,  are  the 
preferred  manner  for  utilization  of  magnitude  regression  results. 
However,  in  this  instance,  a  point  estimate  is  required  by  the  curve 
extension  procedure. 

Thus,  the  maximum  magnitude  derived  from  this  technique  is  : 

Ms  =  6.95, 

Here,  mb(Lg)  has  been  taken  as  equal  to  mb  and  the  Nuttli  and  others 
(1989)  mb  to  Ms  conversion  has  been  usec^ 
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Other  Statistical  Approaches 


These  approaches  make  use  of  extreme-value  theory  (see,  e.g., 
Yegulalp  and  Kuo,  1974  or  Kijko,  1984).  That  theory  assumes  that  the 
occurrence  of  maximum  earthquakes  within  a  given  interval  of  time  is  a 
random  event  and  that  maximum  earthquakes  in  the  future  will  occur  in 
the  same  way  as  those  in  the  past.  In  principle,  this  sounds  ideally  suited 
to  the  task  at  hand.  However,  in  applications  to  date,  it  appears  that  very 
large,  high  quality  data  sets  are  required  for  useful  results  (Coppersmith 
and  others, 1 987).  Knopoff  and  Kagan  (1977)  studied  synthetic  data  sets  to 
show  that  unacceptably  large  errors  could  result  from  extreme-value 
techniques  with  data  bases  similar  to  those  often  encountered  in  practice. 

McGuire  (1977)  investigated  the  use  of  the  sparse  data  sets  available 
for  the  eastern  U.S.  to  estimate  the  maximum  earthquake  by  means  of 
maximum-likelihood  techniques.  He  concluded  that  the  data  were 
inadequate  to  define  with  any  confidence  the  maximum  possible 
earthquake  for  a  given  seismic  zone.  Bender  (1988)  extended  McGuire's 
conclusion  to  state  that,  for  most  real  data  sets,  the  amount  of 
information  available  is  too  small  to  permit  a  reliable  estimate  of 
maximum  magnitude  to  be  obtained,  regardless  of  the  technique  used. 
Accordingly,  this  class  of  estimation  procedures  will  not  be  applied  in 
this  study. 


Applications  of  Fault  Zone  Dimensions 

Magnitude  versus  Fault  Area  Results  for  the  GCSZ 

The  earthquake  foci  within  ±  10  km  of  a  plane  trending  N25°W  are 
assumed  to  be  in  the  GCSZ  and  they  will  be  used  to  estimate  the  spatial 
dimensions  of  the  causal  geologic  fault  zone  structure  (Bollinger  and 
others,  1989b).  The  bases  for  that  assumption  are  :  (1)  The  N25°W  trend  is 
the  approximate  strike  of  the  causal  fault  zone,  (2)  The  dip  of  the  zone  is 
steep,  but  not  necessarily  vertical,  and,  (3)  The  actual  errors  in  the 
hypocentral  locations  may  be  somewhat  larger  than  estimated.  That 
estimate  is  assumed  to  be  for  the  same  geologic  structure  that  was 
involved  in  the  1897  earthquake  sequence,  including  its  Ms  5.9  mainshock, 
as  well  as  the  subsequent  seismic  activity  in  the  Giles  county  locale  up  to 


the  present  time. 

It  is  important  to  note  that,  within  the  past  20  years,  the  'off  zone' 
activity  (Figure  1)  has  included  shocks  of  up  to  a  mb  of  4.6  That  fact  is 
seen  as  being  compatible  with  the  spatial  stress  perturbations  induced 
into  the  volume  surrounding  a  seismic  zone  capable  of  generating  a  Ms  6  or 
larger  earthquake. 

The  areal  extent  of  the  GCSZ  is  described  by  a  small  number  (13)  of 
accurately  located  microearthquake  foci.  They  are  judged  to  be  sufficient 
to  estimate  a  range  for  the  fault  area  by  means  of  a  maximum-minimum 
type  of  approach  that  incorporates  different  assumptions  on  the  shape  of 
the  zone.  Specifically,  six  different  areas  will  be  derived  from  the 
following  : 


*  Two  different  horizontal  lengths,  20  km  and  30  km, 

*  Three  different  vertical  extents,  8  km,  10  km  and  13  km,  and 

*  One  rectangular  shape  and  two  polygonal  shapes. 

The  different  horizontal  and  vertical  dimensions  and  the  different 
configurations  are  necessitated  by  the  inclusion  or  exclusion  of  peripheral 
foci  for  the  purpose  of  estimating  areal  maxima  and  minima  (see  Figure  6 
and  7  ;  Bollinger  and  others,  1989b).  These  6  areas  will  now  be  used  to 
determine  6  magnitudes  whose  mean  value  will  comprise  the  maximum 
magnitude  estimate  for  this  technique. 

Rectangular  Fault  Shape  :  The  vertical  section  (Figure  6)  shows  that 
the  well-constrained  focal  depths  in  the  zone  vary  from  about  5  km  to 
about  15  km.  The  horizontal  extent  is  for  lengths  of  approximately  20  km 
or  30  km  depending  on  whether  or  not  the  two  most  northeasterly  foci,  at 
approximately  15  km  depth,  are  included  or  not.  Assuming  a  simple 
rectangular  shape  from  these  approximate  dimensions  yields  areas  of  200 
sq  km  or  300  sq.  km.  These  areas,  in  turn,  imply  Ms  values  of  6.59  and  6.76 
respectively. 

Polygonal  Fault  Shapes  :  Instead  of  using  the  spatial  distribution  of 
foci  as  a  general  guide  as  was  done  in  the  preceding,  they  can  also  be 
employed  as  the  actual  periphery  of  the  zone.  For  that  senario,  the 
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outermost  hypocenters  are  connected  by  straight  lines  and  the  enclosed 
area  measured.  In  this  instance,  that  procedure  allows  for  three  different 
vertical  dimension  values  that  differ  principally  depending  on  whether  or 
not  the  single  deep  focus  at  18  km  is  included  or  not  (Figure  7).  The 
resultings  areas  are  112,  157,  190,  and  253  sq  km.  The  derived  Ms 
estimates  are  6.34,  6.48,  6.57  and  6.69  respectively. 

The  average  of  the  preceding  six  estimates  is  ; 

Ms  =  6.57, 

and  that  value  serves  as  the  maximum  magnitude  estimate  for  this 
procedure. 

Bonilla  and  others'  (1984)  Magnitude  versus  Fault  Length  relationship 
yields  Ms  =  6.96  for  L  =  20  km  and  Ms  =  7.10  for  L  =  30  km.  It  is  interesting 
to  note  that  these  values  are  closer  to  those  from  the  incremental 
methods  than  the  average  from  the  areal  method. 


Strain  rate  or  Rate  of  Moment  Release  Methods 

These  techniques  have  been  developed  in  seismically  active, 
interplate  regions,  such  as  California,  where  the  active  faults  are 
available  for  geologic  study  by  surface  methods  and  the  strain  rates  are 
high  enough  to  be  measurable  by  geodetic  and  seismic  means.  Such 
conditions  and  data  bases  are  not  available  for  the  eastern  U.S.  and.  thus, 
this  class  of  methods  cannot  be  brought  to  bear  on  the  problem  at  hand. 

Reference  to  a  Global  Data  Base 

The  rationale  here  is  to  substitute  space  for  time  in  an  attempt  to 
overcome  a  short  historical  record  as  in  the  EPRI  study  by  Coppersmith  and 
others  (1987).  Their  results  can  be  employed  as  a  qualitative  tool  to 
assist  in  the  estimation  of  maximum  earthquakes.  Those  results  to  date 
are  ; 

1)  Only  5  great  earthquakes  (M>8)  and  some  20  shocks  larger  than  Ms 
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7  worldwide  in  intraplate  regions  during  historical  time, 

2)  Most  of  those  earthquakes  (68%)  are  at  locations  of  prior 
seismicity. 

3)  Paleozoic  crust  is  far  more  active  crustal  age  province  compared 
to  Precambrian  crust,  and, 

4)  The  horizontal  deviatoric  stress  is  compressive  in  86%  of  the 
cases. 

Unfortunately,  at  its  present  state  of  development,  this  technique  is  not 
suitable  for  application  to  the  problem  at  hand. 

Summary 

Three  estimates  for  the  maximum  magnitude  associated  with  the 
Giles  county,  Virginia,  seismic  zone  have  been  developed.  Those  estimates 
and  the  techniques  employed  to  obtain  them  are  : 

Ms  =  6.9,  from  adding  an  increment  to  the  maximum  historical 
earthquake  in  the  zone, 

Ms  =  6.95,  from  extension  of  the  magnitude  recurrence  curve  for 
the  zone,  and 

Ms  =  6.57,  from  estimates  of  the  area  of  the  zone. 

These  values  are  to  be  given  equal  weight  and  can  be  employed  in  a  number 
of  ways.  If  a  single  value  is  required,  then  the  average,  rounded  to  the 
nearest  one-tenth,  should  be  used.  That  value  is  : 

Ms  =  6.8  . 

The  mb  equivalent  is  6.3.  If  multiple  values  can  be  accommodated  the  two 
extreme  values  can  be  utilized. 
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APPENDIX  A. 


Earthquake  Catalog  for  the  Giles  County,  Virginia,  Seismic  Zone  - 
All  earthquakes  :  1 876  through  1 988 
within 

A  circle  of  radius  50  km  centered  at  37.25°  -  80.75° 
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13-FEB-89  14: 19:30  for  the  prog-am  RWGEN 

Input  data  from  file  .  GCO.  GEN 
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9 

N 

81.  1 

W 

18790901 

12 

36. 

9 

N 

81.  1 

W 

18850202 

1210 

36. 

9 

N 

81.  1 

W 

18970503 
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1 

N 
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W 

18970503 

19 
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N 
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W 
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N 
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W 

18970503 

23 

3-' 

1 

N 

80.  7 

w 

18970531 

1858 

37 

3 

N 

80.  7 

w 

18970629 

03 

37. 

3 

N 

80.  7 

w 

18970904 

1 1 

36. 

9 

N 

81.  1 

w 

18971022 

0320 

36. 

9 

N 

81.  1 

w 

18980205 

20 

37. 

0 

N 

81.  0 

w 

18980206 

02 

37. 

0 

N 

81.  0 

w 

1B981125 

20 

37 

0 

N 

81. 0 

w 

18990213 

0930 

37 

0 

N 

81. 0 

w 

19020518 

04 

37 

3 

N 

80.  6 

w 

19170419 

37. 

0 

N 

81.  0 

w 

19590423  205839.  5  37. 395N  80.  682W  1.  0 

19590707  2317  37  3  N  80.  7  W 

19590821  1720  37  3  N  80.  7  W 

19680308  053815.7  37.  281N  80.  774W  8.0 

19691120  010009.3  37  449N  80. 932W  30 

19740530  212835.  3  37.  457N  80.  540W  5.  0 

19750307  124513.  5  37  32  N  80.  48  W  5.  0 

19751111  081037.6  37,  217M  80.  892W  1.0 

19760703  205345.  8  37.  32  N  81 .  13  W  1.0 

19780128  231323.  4  37.  22eN  80.  747W  4.  5 

19780510  041910.9  37.  294N  80.  729 W  17.6 
19780525  083025.  I  37. OOON  80.  794W  12.  1 
19780728  083940.7  37  337N  80.  690W  11.8 
19780830  021938  2  37  362N  80  668W  8.  4 

19800218  035855. 2  37.  42SN  80  593W  14  6 

19800410  223315  5  37  496N  81  lllW  0.  2 

19801202  074738.  2  37.  414N  80  546W  13,  1 

19810824  115011.4  36.  949N  80.  742W  16.1 
19811112  062414  1  37  235N  80.  744W  7.9 
19811204  023556.5  36.  999N  80.  754W  5,9 
19820518  031633  8  37. 130N  80. 497W  11.0 
19830108  155355.9  37  331N  80. 614W  67 

19830125  203858.  0  37,  386N  80.  509W  16.  8 
19830420  180956.  4  37. 356N  80.  825W  11.5 
19830517  020247.  6  37  254N  80. 733W  7.  0 

19830526  010444. 9  37  507N  80. 315W  8.  6 

19830710  140539.5  37. 272N  80.  752W  7.7 

19831113  165106.6  37.  566N  80.  759W  11.8 
19831113  175049.8  37  572N  80.  760W  14  2 
19831209  001158.0  37  200N  80.  786W  12.4 
19840311  040139.0  37.  470N  80.  890W  1.  1 

19840702  195138.  6  37  278N  80.  725W  10.  8 

198411)7  031728.3  37  266N  80.  727W  10.3 
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There  have  been  50  events  listed  so  far. 
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1  Date 

lOT  (UCT) : 

Hypocenter 

iError 

(km)  ISr! 

Magnitudes 

1  Int 

Lab 

YearMoDy 

HrMnSec 

Lat 

Lon 

Depth  ERH 

ERZ 

Magi 

Mag2  Mag3 

S 

VA 

19850610 

122238. 3 

37. 248N 

80.  48SW 

11.1  0.  9 

2.  0  IV 

3.  2NV 

2.  8DV 

4V 

UIV 

19850614 

075710. 2 

37.  534N 

81.  020W 

2.  4  3.  6 

4.  4  IV 

0.  0 

0.  SDV 

VA 

19860326 

163623.  9 

37. 245N 

00  494W 

11.9  1.1 

2.  4  IV 

2.  9 

2.  9DV 

4V 

There  are  53  events  in  this  listing. 


SOURCE  CODES: 

B  -  Bollinger.  1975.  Southeastern  U.  S.  Catalog  1754-1974. 

E  ~  Earth  Physics  Branch.  Canadian  catalog. 

G  -  USGS  -  State  Seismicity  Maps  (Stover /Reagor  et  al.  ). 

I  -  EPRI  Catalog  (8  July  1906). 

N  -  Neilsen,  1982  (Stanford  Data  Base. . .  ). 

R  -  Barstou)  et  al.  .  1981  (Rondout  Asso.  ),  NUREG/CR- 1 577, 

S  -  Street  and  Turcotte.  1977.  BSSA,  67.  pp.  599-614, 

T  -  Reinbold  and  Johnston  (TEIC),  1986,  USGS  Final  Ropt.  , 

U  -  Earthquake  History  of  the  U.  S.  /U.  S.  Earthquakes. 

-  SEUSSN  Bulletins  (Va.  Tech  Publication), 

Y  -  Felt  area  only;  value  is  the  average  of  those  found  in  G  and  R  above, 

Z  -  Felt  area  only;  value  is  the  average  of  those  found  in  U  and  R  above. 

LOCATION  CODES: 

H  -  Historical  Location  (from  intens i ty / f e 1 t  area  data), 

I  -  Instrumental  Location, 

MAGNITUDE  CODES: 

B  -  mb  from  Baysian  estimate  (Veneziano  it  VanDyck,  1984), 

C  -  mb  from  intensity  and  felt  area  (Sibol  et  al. ,  1987) 

D  -  Md  from  duration  or  coda  length, 

F  -  mb  from  felt  area/attenuation  data, 

I  -  mb  from  intensity  data 

L  -  ML  (Richter,  1958), 

M  -  mb  determined  from  modified  i ns tr umen t s / f or mu  1 1 , 

N  -  mb  from  Lg  wave  data  (Nuttli,  1V73). 

0  -  m3Hz  (Lawson,  et  al.  ,  1979  -  Oklahoma  ear th qua k es ) , 

P  -  mb  from  P  wave  data  (Gutenberg  and  Richter,  1956), 

S  -  MS  (Bath.  1966;  Gutenberg,  1945), 

X  -  Magnitude  of  unknown  type. 


E28 


APPENDIX  B. 


Earthquake  Catalog  for  the  Giles  County,  Virginia,  Seismic  Zone  - 
All  earthquakes  whose  hypocentral  error  estimates  are  <  ±10  km 

and 

whose  epicenters  are  within  a  rectangle  oriented  N25°E 
with  dimensions  of  20  km  by  41  km. 


E29 


13-FEB-89  14  19: 52  for  the  program  RWGEN 

Input  data  from  file .  CCOSMALL.  CEN 


Page 


Lab 

!  Date 
YearMoDy 

:  QT  (UCT) : 
HrMnSec 

Hypocenter 

Lat  Lon  Depth 

1  Error 
ERH 

(  km ) 
LRZ 

ISr  ! 

Magn i tud 
Magi  Mag2 

VA 

19780128 

231323. 4 

37.  228N 

80. 747W 

4.  5 

5.  9 

3.  0 

IV 

1.  6 

1.  6DV 

VA 

19780510 

041910. 9 

37.  294N 

80.  729W 

17.  6 

2.  7 

2.  4 

IV 

0.  a 

0.  8DV 

VA 

19780728 

083940. 7 

37.  337N 

80.  690W 

11.8 

4  9 

8.  1 

IV 

0.  6 

0  6DV 

VA 

19780830 

021938. 2 

37  362N 

80. 668U 

8.  4 

3.  1 

6.  4 

IV 

0.  5 

0.  5DV 

VA 

19800218 

035855.  2 

37. 428N 

80. 593W 

14.  6 

0.  6 

0.  7 

IV 

0.  6 

0.  6DV 

VA 

19801202 

074738.  2 

37.  414N 

80.  546W 

13.  1 

1.  1 

2.  3 

IV 

0.  5 

0.  5DV 

VA 

19811112 

062414. 1 

37. 235N 

80. 744 W 

7.  9 

2.  6 

8.  0 

IV 

0.  8 

0.  8DV 

VA 

19830108 

155355.  9 

37. 331N 

80.  614W 

6.  7 

1.  3 

5.  7 

IV 

1.  3 

1.  3DV 

VA 

19830517 

020247.  6 

37. 254N 

80.  733 W 

7.  0 

0.  9 

2.  9 

IV 

0.  0 

0.  ODV 

VA 

19830710 

140539.  5 

37.  272N 

80.  752W 

7.  7 

0  7 

1.  9 

IV 

1.  1 

1.  IDV 

VA 

19831209 

001158. 0 

37. 200N 

80.  78614 

12.  4 

0.  5 

0.  a 

IV 

1.  5 

1.  5DV 

VA 

19840702 

195138. 6 

37. 278N 

80. 725W 

10.  8 

0.  5 

0.  9 

IV 

1.  5 

1.  5DV 

VA 

19841117 

031728. 3 

37.  266N 

80.  727W 

10.  3 

0.  4 

1.  0 

IV 

0.  0 

0.  ODV 

Mag  3 


!  Int 
S 


There  are 


13  events  in  this  listing. 


SOURCE  CODES: 

B  -  Bollingeri  1975i  Southeastern  U.  S.  Catalog  1754-1974. 

E  -  Earth  Phgsics  Branch.  Canadian  catalog. 

Q  -  uses  -  State  Seismicity  Maps  (Stover /Reagor  ct  al.  >. 

I  -  EPRI  Catalog  <8  July  1986), 

N  -  Neilsen.  1982  (Stanford  Data  Base.  .  .  ). 

R  -  Barstow  et  al.  .  1981  (Rondout  Asso.  ).  NUREO/CR-1 577. 

S  -  Street  and  Turcotte.  1977.  8SSA,  67.  pp.  599-614. 

T  -  Reinbold  and  Johnston  (TEIC),  1986.  USGS  Final  Rept.  . 

U  -  Earthquake  History  of  the  U.  S.  /U.  S.  Earthquakes. 

V  -  SEUSSN  Bulletins  (Va.  Tech  Publication). 

Y  -  Felt  area  only;  value  is  the  average  of  those  found  in  C  and  R  above. 

Z  -  Felt  area  only;  value  is  the  average  of  those  found  in  U  and  R  above. 

LOCATION  CODES; 

H  -  Historical  Location  (from  intens i ty /f e 1 t  area  data). 

I  -  Instrumental  Location. 

MAGNITUDE  CODES. 

B  -  mb  from  Baysian  estimate  (Veneziano  Sc  VanDyck.  1984). 

C  -  mb  from  intensity  and  felt  area  (Sibol  et  al. .  1987) 

D  -  Md  from  duration  or  coda  length. 

F  -  mb  from  felt  arca/attenuation  data, 

I  -  mb  from  intensity  data 

L  -  ML  (Richter,  1958), 

M  -  mb  determined  from  modified  instruments/f ormol i. 

N  -  mb  from  Lg  wave  data  (Nuttli,  1973), 

0  -  m3Hz  (Lawson,  et  al.  ,  1979  -  Oklahoma  ear thquakes ) , 

P  -  mb  from  P  wave  data  (Gutenberg  and  Richter.  1956), 

S  -  MS  (Bath.  1966;  Gutenberg,  1945), 

X  -  Magnitude  of  unknown  type. 


APPENDIX  F: 

RECOMMENDED  ACCELEROGRAMS  AND  RESPONSE  SPECTRA 

From  California  Institute  of  Technology, 
Strong  Motion  Earthquake  Catalogue,  1971  to  1975 


FI 


SRN  FERNRNDO  ERRTHQURKE  FEB  9.  1971  -  0600  PST 

10198  71.069.0  GRIFFITH  PRRK  OBSERVRTORY.  MOON  ROOM.  LOS  RNGELES.  CRL .  COMP  SOOW 

o  PERK  VRl  UES  :  RCCEL  =  -176.9  CM/SEC/SEC  VELOCITY  =  -20.2  CM/SEC  DI5PL  -  7.3  CM 


o  o 

lD  fO 
CM  1 


O  OO 
ro  ' 


33S/33S/W3 

N0I18y333338 


33S/W3 

AiI3013A 


W3 

iN3W33U1dSia 


SECONDS 


OIII 


F3 


VELOCITY  (m./sec) 


RESPONSE  SPECTRUM 


SP'J  FER'iS'.QJ  Er.RTH^L'fiaS  FCB  9.  197;  -  0600  P5T 
niri93  71.C39.0  (ifiifFiTH  pr,;r  cBSfovaTCFtT,  n.-ch  lcs  p'-'-c.cc.  cfl.  crf  s:rH 
DP.visG  vrlues  0.  ?.  s,  to  f*.o  ro  peaces:  cfiiTicsi. 


PERIOD  (secs) 


RELfiTIVE  VELOCITY  RESPONSE  SPECTRUM 

SflN  FERNANDO  EARTHQUAKE  FEB  9.  1971  -  0600  PST 
11IG106  71.018.0  CALTECH  SEISMOLOGICAL  LAB.,  PASADENA.  CAL.  COMP  S90W 

DAMPING  VALUES  ARE  0.  2.  5.  10  AND  20  PERCENT  OF  CRITICAL 


PERIOD  -  SECONDS 


RESPONSE  SPECTRUM 


SfIN  FERNANDO  ERRTHQURKE  FEB  9.  1971  -  0600  PST 

IIIGtOe  71.018.0  CALTECH  SEISHOLOGICfiL  LRB.,  PRSROENfl,  CflL.  COHP  S90H 
DRMPING  VRLUES  RRE  0.  2.  5.  10  RND  20  PERCENT  OF  CRITJCRL 


PERIOD  (secs) 


F7 


SflN  FERNANDO  EARTHQUAKE  EEB  9.  1971  -  0600  P5T 
1IP221  71. 150. Q  SANTA  ANITA  RESERVOIR,  ARCADIA.  LAL .  COMP  NOSE 

©PEAK  VALUES  :  ACCEL  =  -137.7  CM/SEC/SEC  VELOCITY  =  S.O  CM/SEC  DI5PL  = 


SECONDS 


RELATIVE  VELOCITY  RESPONSE  SPECTRUM 


F9 


FtniGD  -  SECONDS 


VELOCITY  (in. /sec) 


RESPONSE  SPECTRUM 
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